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EXECUTIVE  SUMMARY 


The  damaging  effects  from  an  uncontained  aircraft  turbine  engine  failure  can  be  catastrophic.  As 
a  result,  the  Federal  Aviation  Administration  (FAA)  has  commissioned  a  program  to  mitigate  the 
damaging  effects  of  such  an  event.  The  Uncontained  Engine  Debris  Mitigation  Program  will 
involve  both  industry  and  government  to  determine  possible  engineering  solutions  to  this 
problem. 

As  part  of  this  program,  the  Naval  Air  Warfare  Center  Weapons  Division  (NAWCWPNS)  has 
been  tasked  to  evaluate  ballistic  damage  analysis  tools  and  techniques  that  are  currently  in  use  by 
the  defense  community.  The  intent  is  to  determine  their  applicability  in  predicting  the  damaging 
effects  from  an  uncontained  engine  failure.  This  report  documents  testing  that  was  conducted 
and  the  evaluation  of  several  empirical  penetration  equations  under  the  circumstances  present 
during  engine  failure  events. 

The  data  generated  under  this  effort  showed  that  variations  in  the  velocity,  orientation,  and  shape 
of  the  debris  results  in  differing  failure  modes  of  the  targets.  When  plugging  failures  of  the 
targets  occur,  the  results  of  the  defense  equations  are  quite  good.  When  petaling  failures  occur 
the  prediction  accuracy  was  degraded.  This  resulted  in  the  development  of  an  FAA  Energy 
Equation  through  numerical  curvefitting  of  the  test  data  for  2024  T3  aluminum  targets. 

Additional  testing  is  planned  to  evaluate  the  penetration  equations  performance  with  real  aircraft 
structure  and  varying  obliquity  angles  at  impact  under  phase  II  testing. 
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1.  INTRODUCTION. 


1.1  BACKGROUND  AND  OVERVIEW. 

An  uncontained  engine  failure  precipitated  the  loss  of  control  of  a  DC- 10  at  cruise  altitude  over 
the  Midwestern  United  States  in  1989.  The  aircraft  crashed  while  attempting  to  land  at  Sioux 
City  Airport.  The  accident  investigation  revealed  that  multiple  fragments  from  the  engine  failure 
damaged  all  of  the  redundant  flight  control  systems  resulting  in  a  loss  of  the  hydraulic  systems. 
The  aircrew  maintained  a  modicum  of  control  by  varying  the  power  level  on  the  remaining  wing 
mounted  engines  saving  many  lives. 

The  National  Transportation  Safety  Board  (NTSB)  produced  several  recommendations  as  a  result 
of  that  accident  and  this  effort  continues  to  support  the  Federal  Aviation  Administration’s  (FAA) 
response  to  those  recommendations  which  are  intended  to  develop  revised  certification 
requirements  to  minimize  potential  for  catastrophic  failures  from  multiple  fragments  during  an 
uncontained  engine  failure.  This  effort  is  being  performed  under  the  advice  of  the  Aviation 
Rulemaking  Advisory  Committee  for  Power  Plant  Installation  and  Harmonization  Working 
Group  with  industry  and  government  to  develop  analytical  methods  to  evaluate  the  problem  and 
mitigate  the  consequences. 

As  part  of  this  program,  the  Naval  Air  Warfare  Center  Weapons  Division  (NAWCWPNS)  has 
been  tasked  to  evaluate  ballistic  damage  analysis  tools  and  techniques  which  are  currently  in  use 
by  the  defense  community.  The  intent  is  to  determine  their  applicability  in  predicting  the 
damaging  effects  from  an  uncontained  engine  failure.  This  report  documents  testing  which  was 
conducted  in  the  evaluation  of  these  equations 


1.2  THE  JTCG/ME  PENETRATION  EQUATIONS. 

Among  the  analysis  tools  currently  in  use  are  the  Joint  Technical  Coordinating  Group  for 
Munitions  Effectiveness  (JTCG/ME)  penetration  equations  [1].  These  equations  were  originally 
developed  to  predict  the  penetration  and  residual  velocity  of  ballistic  weapon  projectiles  upon 
impact  with  various  materials.  Of  particular  interest  are  the  ballistic  prediction  equations  which 
were  developed  for  warhead  fragments,  as  these  fragments  are  closest  in  geometric  description  to 
jet  engine  blade  fragments.  The  JTCG/ME  penetration  equations  are  listed  as  follows. 
Derivation  of  these  equations  is  covered  in  appendix  A. 
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=  debris  initial  velocity 

=  debris  ballistic  limit  velocity  (Penetration  occurs  50%  of  the  time.) 
=  debris  residual  velocity 
=  debris  specific  weight  (debris  weight/volume) 

=  plate  specific  weight  (plate  weight/volume) 

=  plate  thickness 

=  debris  presented  area  along  direction  of  travel 

=  debris  obliquity  angle  at  impact 
=  debris  weight 
=  100  grains 


Empirical  Constants:  Cbf,bf,h,f 


The  penetration  equations  have  been  incorporated  into  a  computer  model  as  one  possible  analysis 
tool  for  aircraft  designers.  The  Uncontained  Debris  Model  (UDM)  focuses  on  blade  fragment 
debris  as  this  type  of  debris  is  most  likely  to  result  from  an  uncontained  engine  event.  According 
to  extensive  investigation,  including  the  development  of  an  extensive  uncontained  failure 
database,  these  are  the  most  likely  particles  which  would  result  in  multiple  damage  sites  to  an 
aircraft.  The  debris  model  identifies  the  debris  type  by  component  (fan  blade,  turbine  blade, 
compressor  blade)  and  defines  the  most  likely  size,  mass,  and  velocity  for  each  debris  fragment. 

2.  TEST  OBJECTIVES. 

The  objective  of  these  tests  was  to  collect  fragment  penetration  data  to  validate  and/or  refine  the 
penetration  equations  used  to  develop  the  Uncontained  Debris  Model.  Testing  was  conducted 
using  simulated  and  actual  engine  debris  fragments  against  various  thicknesses  of  aluminum 
plates  and  aircraft  skin  structures. 


3.  APPROACH. 

3.1  TEST  OVERVIEW. 

The  data  required  for  validation  of  the  model  would  be  penetration  data  resulting  from  the  impact 
of  engine  blade  fragments  against  representative  aircraft  fuselage  material.  This  data  would 
include  the  impact  and  residual  velocities  of  the  fragments,  along  with  the  impact  orientation  of 
each  fragment  with  respect  to  the  target. 

To  achieve  the  collection  of  this  data,  it  would  be  necessary  to  propel  representative  engine  blade 
fragments  at  velocities  representative  of  an  uncontained  engine  event  (approximately  200- 
800  ft/s).  It  was  decided  to  utilize  the  Missile  Impact  Kinetic  Energy  Simulator  (MIKES)  airgun 
to  propel  the  engine  blade  fragments. 

The  MIKES  gun  is  a  recent  addition  to  the  NAWCWPNS  Weapons  Survivability  Laboratory 
(WSL).  Although  originally  developed  as  a  means  to  propel  man  portable  missiles,  its  large 
barrel  diameter  (6.07")  would  be  capable  of  accommodating  the  larger  representative  fan  blade 
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fragments.  Figure  1  illustrates  the  test  site  setup.  Specifics  of  this  setup  are  outlined  in  the 
detailed  test  plan  included  in  appendix  D. 


FIGURE  1.  TEST  SITE  SETUP  WITH  MIKES  GUN 


The  fragments  were  propelled  from  the  gun  by  means  of  a  sabot.  At  the  end  of  the  gun  barrel,  a 
sabot  stopper  was  placed  to  stop  the  sabot  but  allow  the  fragment  to  continue  travel  unimpeded. 
The  original  design  as  outlined  in  the  test  plan  for  the  sabot  and  sabot  stopper  proved  impractical. 
The  sabot  would  shatter  within  the  stopper  allowing  sabot  fragments  to  travel  along  with  the 
engine  fragment. 


Alternative  designs  for  the  sabot  and  stopper  were  developed  which  utilized  a  rubber  sabot  along 
with  an  energy  absorbing  sabot  stopper.  Details  of  these  designs  are  illustrated  in  figures  2  and 
3.  The  new  designs  were  successful  in  stopping  the  sabot  and  allowing  unimpeded  travel  of  the 
engine  fragment.  The  rubber  sabot  additionally  proved  durable,  allowing  its  partial  reuse  for 
successive  shots. 


FIGURE  2.  RUBBER  SABOT 


Rubber  Bumpers 


FIGURE  3.  SABOT  CATCHER 


3.2  TEST  SEQUENCE. 

The  test  matrix  outlined  in  the  test  plan  called  for  nine  different  test  series  to  be  conducted. 
These  series  propelled  various  types  and  sizes  of  engine  blade  fragments  at  current  and  proposed 
fuselage  materials.  A  total  of  eighty-four  shots  were  conducted  during  this  testing.  Six  of  these 
shots  failed  to  produce  data,  mainly  due  to  target  misses.  For  those  shots  which  were  successful, 
impact  and  residual  velocities  of  the  fragments  were  measured  through  the  use  of  high-speed  film 
and  breakpapers.  Impact  orientation  of  the  fragment  was  also  determined  from  high-speed  film 
views  within  two  perpendicular  planes. 

Fragments  utilized  included  engine  fan,  turbine,  and  compressor  blade  fragments.  Figure  4 
illustrates  the  various  fragments  utilized  in  this  testing.  An  additional  fragment  was  utilized 
which  was  similar  to  that  used  during  the  1997  SRI  International  Advanced  Armor  Technology 
impact  experiments  [2],  All  fragments  were  solid  titanium  with  the  exception  of  the  turbine 
blade  fragments  which  were  composed  of  steel. 

3.2.1  Test  Series  #1. 

This  test  series  verified  the  effectiveness  of  the  sabots  in  achieving  end  on  and  tumbled  impacts 
with  the  target.  The  accuracy  of  fragment  velocity  calculations  were  also  verified  by  comparison 
of  breakpaper  computations  with  high-speed  camera  visual  calculations.  Finally,  proper 
operation  of  all  test  systems  were  verified  during  these  tests. 

During  the  conduct  of  this  series,  it  immediately  became  apparent  that  achieving  impact 
orientations  as  specified  for  each  test  series  would  be  difficult.  The  engine  blade  fragments 
tumbled  severely  soon  after  leaving  the  muzzle  of  the  gun.  It  was  determined  from  high-speed 
camera  footage  that  this  was  primarily  due  to  the  airfoil  shape  of  these  fragments.  High-speed 
travel  through  the  air  created  a  large  lift  vector  which  was  off-center  from  the  fragment’s  center 
of  gravity.  This  moment  resulted  in  the  fragment’s  tumbling. 
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As  this  tumbling  would  also  likely  occur  during  an  engine  burst  event,  it  was  decided  to  proceed 
with  a  change  to  the  test  matrix.  All  fragments  would  be  propelled  at  an  initial  end-on 
orientation  within  the  sabot.  Fragment  tumbling  and  the  resulting  impact  orientation  would  be 
observed  and  measured  from  the  high-speed  camera  footage.  Analysis  of  the  penetration 
equation’s  prediction  effectiveness  would  be  based  on  each  fragment’s  orientation  at  the  moment 
of  impact. 

The  series  also  verified  the  effectiveness  of  the  velocity  measurement  methods  proposed  in  the 
test  plan.  It  was  noted  that  fragment  tumble  often  created  a  deflection  of  the  fragment  along  the 
intended  shotline.  Penetration  of  the  fragment  through  the  target  created  additional  deflections. 
As  such,  the  fragments  sometimes  missed  the  target  breakpaper.  This  breakpaper  was  also  easily 
damaged  by  debris  in  the  gun’s  air  blast  prior  to  fragment  penetration.  Additionally,  the  residual 
velocity  breakpapers  placed  well  behind  the  target  were  rarely  impacted  due  to  postimpact 
deflection  of  the  fragments.  For  shots  which  did  achieve  proper  breakpaper  operation, 
breakpaper  velocity  calculations  were  compared  with  velocity  calculations  from  the  high-speed 
film  footage.  Both  were  in  close  agreement. 

Due  to  the  difficulties  encountered  with  use  of  the  breakpapers,  it  was  decided  to  utilize  the 
target  breakpaper  as  a  means  to  quickly  verify  that  the  desired  impact  velocity  range  was 
achieved  for  each  test  shot.  As  a  cost  saving  measure,  the  residual  velocity  breakpapers  were  not 
utilized  for  subsequent  tests.  It  was  also  decided  that  the  main  source  of  velocity  data  for  these 
tests  would  be  the  high-speed  camera  footage  due  to  its  close  agreement  with  the  breakpapers. 

3.2.2  Test  Series  #2. 

Test  series  #2  fired  simulated  blade  fragments  similar  to  those  utilized  during  the  SRI 
International  Advanced  Armor  Technology  impact  experiments.  The  dimensions  of  these 
rectangular  titanium  fragments  were  1.4"  length  x  1.0"  width  x  0.25"  thickness.  Each  fragment 
weighed  approximately  25  grams.  The  impact  target  was  0.05"  2024-T3  aluminum  target  plates. 
This  test  series  provided  baseline  impact  resistance  data  for  2024-T3  aluminum.  The  series  also 
allowed  comparison  with  the  accumulated  data  from  the  SRI  experiments.  A  total  of  eleven 
shots  were  conducted  for  this  test  series.  However,  only  five  of  these  tests  provided  useful  data. 

Various  difficulties  were  encountered  during  this  test  series  which  degraded  the  acquisition  of 
data.  These  included  three  firings  which  missed  the  target  at  the  10-foot  standoff  distance 
prescribed  in  the  test  plan.  As  a  result,  the  target  was  moved  to  a  new  standoff  of  5  feet. 

One  firing  was  conducted  with  overexposed  high-speed  film.  Also  the  small  simulated  fragment 
was  indiscernible  in  the  high-speed  coverage  of  two  additional  test  firings.  This  was  due  to  the 
fragment’s  speed  and  orientation  with  respect  to  the  camera  for  these  firings.  As  a  result, 
fragment  velocities  and  impact  orientation  were  not  accurately  determined  for  these  shots.  As 
the  model  predictions  are  highly  dependent  on  these  values,  these  three  firings  did  not  provide 
useful  data. 
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3.2.3  Test  Series  #3  Through  7. 


These  test  series  impacted  actual  engine  blade  fragments  against  both  0.05"  and  0.07"  2024-T3 
aluminum  target  plates.  Characteristics  of  these  fragments  (dimensions,  material,  weight,  and 
impact  velocity/orientation)  were  those  which  were  most  likely  to  be  experienced  during  an 
engine  burst.  Test  series  3  and  4  utilized  representative  titanium  fan  blade  fragments  against 
0.05"  and  0.07"  thickness  2024  T3  aluminum  targets  respectively.  Test  series  5  and  7  examined 
titanium  compressor  blade  impacts  while  test  series  6  utilized  steel  turbine  blade  fragments.  Test 
series  5,  6,  and  7  all  utilized  0.05"  2024  T3  aluminum  as  targets  for  the  fragments. 

3.2.4  Test  Series  #8. 

This  test  series  examined  the  reaction  of  engine  fan  blade  impacts  against  an  actual  engine 
cowling.  The  cowlings  were  composed  of  0.05"  7075  T6  aluminum.  Comparison  was  made 
with  the  data  accumulated  from  the  previously  conducted  tests.  This  comparison  would 
determine  consistency  of  the  data  between  the  use  of  new  sheet  metal  targets  and  aged  engine 
cowling  material. 

3.2.5  Test  Series  #9. 

This  final  test  series  was  to  examine  fragment  impacts  against  Kevlar  fabric  target  material.  Due 
to  difficulty  in  acquiring  a  sufficient  quantity  of  Kevlar  material,  Zylon  material  was  substituted 
in  its  place.  A  total  of  twelve  test  shots  were  conducted  against  the  Zylon  targets  using  each  of 
the  represented  fragments.  The  data  from  these  tests  will  be  included  in  a  separate  report  from 
SRI  International. 

4.  DATA  ANALYSIS. 

To  evaluate  the  prediction  accuracy  of  the  JTCG/ME  Penetration  equations,  it  would  be 
necessary  to  accurately  determine  the  velocity  of  the  fragment  at  impact  along  with  the 
fragment’s  postimpact  residual  velocity.  Additionally,  the  presented  area  of  the  fragment  at 
impact  would  need  to  be  determined  as  input  into  the  penetration  equations.  High-speed  film 
camera  coverage  of  each  impact  was  utilized  to  record  this  data. 

A  high-speed  film  record  of  each  shot  was  acquired  through  the  use  of  two  high-speed  film 
cameras.  These  cameras  were  placed  to  provide  views  of  the  impact  in  two  perpendicular  planes. 
These  views  included  perpendicular  side  and  top  views  of  the  target. 

4.1  VELOCITY  CALCULATION. 

Through  the  use  of  a  film  motion  analyzer,  it  is  possible  to  accurately  determine  the  film’s  frame 
rate  at  the  moment  of  impact.  The  analyzer  also  has  movable  crosshairs  and  provides  Cartesian 
coordinate  readouts  for  the  center  of  the  crosshairs.  Film  frames  can  be  incrementally  stepped 
and  viewed. 


7 


To  determine  fragment  velocities,  the  analyzer  was  utilized  to  determine  the  distance  traveled  by 
the  fragment  along  with  the  number  of  film  frames  required  to  traverse  this  distance.  The 
velocity  could  then  be  calculated  by  the  equation: 

Velocity  =  Camera  Frame  Rate  x  Distance  Traveled/#  of  Frames 

The  frame  rate  for  the  high-speed  cameras  were  preset  to  2,000  frames  per  second  for  all  events. 
As  the  film  cannot  instantaneously  accelerate  to  this  frame  rate,  there  was  some  error  between  the 
desired  frame  rate  and  the  actual  film  frame  rate  at  the  moment  of  the  event.  As  an  aid  in 
determining  the  actual  film  frame  rate,  timing  markers  are  placed  on  the  film  by  the  camera  at  1- 
millisecond  intervals.  The  actual  camera  frame  rate  is  determined  by  counting  both  the  number 
of  frames  and  timing  markers  near  the  moment  of  the  event.  The  frame  rate  is  then  calculated  by 
dividing  the  number  of  frames  by  the  number  of  1-ms  timing  markers. 

The  distance  traveled  by  the  fragment  during  the  event  is  determined  through  the  use  of  the 
analyzer’s  crosshairs.  It  is  first  necessary  to  determine  a  scaling  factor  which  equates  the 
analyzer’s  crosshair  coordinate  readout  with  the  actual  distances  viewed  on  the  film.  This  is 
accomplished  by  equating  the  readout  to  visual  reference  distances  viewed  on  the  film.  For  this 
purpose,  scaled  gridboards  with  one-foot  reference  stripes  were  placed  beside  and  below  the 
target. 

As  most  of  the  fragments  tumbled  during  the  test,  the  crosshairs  are  placed  on  a  visual  estimate 
of  the  fragment  center  of  gravity.  A  Cartesian  reading  is  taken  of  the  fragment  a  few  frames  prior 
to  impact.  A  second  reading  is  taken  one  frame  prior  to  impact.  The  traversed  distance  is  then 
calculated  by  multiplying  the  scaling  factor  by  the  difference  of  the  first  and  second  reading. 

4.2  ORIENTATION. 

4.2.1  Presented  Area. 

The  presented  area  of  the  fragment  to  the  target  is  a  critical  input  value  to  the  penetration 
equations.  An  accurate  determination  of  this  value  is  necessary  to  confidently  evaluate  the 
predictions  provided  by  the  equations.  The  film  motion  analyzer  was  also  utilized  to  determine 
this  value. 

The  high-speed  cameras  were  placed  to  provide  two  views  of  the  fragment  in  two  perpendicular 
planes.  These  views  included  a  side  view  and  a  top  view  of  the  fragment  in  the  vicinity  of  the 
target.  The  film  frame  utilized  in  calculating  the  presented  area  for  each  test  was  the  immediate 
frame  just  prior  to  fragment  impact. 

These  two  views  allowed  calculation  of  the  fragment  presented  area  in  the  third  plane  of  the 
target.  This  was  accomplished  by  assuming  that  the  fragment  was  rectangular  in  shape. 
Cartesian  coordinates  for  the  four  comers  of  the  fragment  were  read  from  the  film  analyzer. 
These  coordinates  were  then  converted  into  width  and  length  vectors  for  the  fragment.  The 
presented  area  of  the  fragment  in  all  three  (top,  side,  and  frontal)  planes  could  then  be  calculated 
by  the  cross  product  of  these  two  vectors. 
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As  the  top  view  had  inherent  parallax  errors  due  to  its  close  proximity  to  the  target, 
determination  of  the  length  and  width  vectors  utilizing  visual  scaling  factors  would  have  incurred 
additional  errors.  Instead,  vector  angles  were  used  to  calculate  the  width  and  length  vectors  as 
described  in  figure  5. 


e  =( 

w,t  \+ 


ArcTan  ( (Xb,s-Xd,s)  /  (Zb,s-Zd,s) ) 
f  ArcTan  ( (Xa,s-Xc,s)  /  (Za,s-Zc,s) ) 

_  /  ArcTan  ( (Xc,s-Xd,s)  /  (Zc,s-Zd,s) ) 
ArcTan  ( (Xa,s-Xb,s)  /  (Za,s-Zb,s) ) 


TOP  VIEW 


)  /  2 

)/. 


Length  Vector:  L  =  Length  ( sin  0Lt  i  +  sin  0Ls  j  +  sin  0Ls  k) 

Width  Vector:  W  =  Width  (sin  ©W(  i  +  sin  0Ws  j  +  sin  ©Ws  k) 

Area  Vector:  A  =  L  X  (Cross)  W 

Presented  Area:  Ap  =  k  component  of  Area  Vector  A 

or 

Ap  =  Length  X  Width  (Sin  ©uSin  0w,s  Sin  9§in 


FIGURE  5.  CALCULATION  OF  FRAGMENT  PRESENTED  AREA 

It  should  be  noted  that  the  calculated  value  of  presented  area  may  be  negative,  depending  on  the 
fragment’s  orientation  at  impact.  If  negative,  the  absolute  value  of  this  term  is  utilized  as  input 
to  the  penetration  equations. 


9 


Figure  5  illustrates  the  calculation  of  the  presented  area  based  on  the  fragment’s  side  dimensions. 
The  edges  of  the  fragment’s  perimeter  may  also  contribute  to  presented  area,  especially  at  near 
edge-on  impacts.  Calculation  of  the  edge’s  contributions  to  presented  area  is  a  bit  more  involved 
than  displayed  in  figure  5.  These  areas  may  be  calculated  by  assuming  that  the  rectangle  is  of  a 
specified  thickness.  A  thickness  vector  T  may  be  calculated  by  assuming  perpendicularity  to  the 
length  and  width  vectors  L  and  W.  As  a  definition  of  perpendicularity,  the  dot  product  of  L  and 
T,  and  W  and  T  will  be  zero.  By  applying  this  assumption  with  much  math  (which  is  not 
included  here),  the  resulting  thickness  vector  may  be  calculated  by  the  following  relations: 


P  = 


(wk.it)  /( h 

V  Wj  Lj  J  J  \  Lj 


Tk 


2 


T  j  —  P  x  Tk 

Ti=-^jxLj+TkxLkyLi 
T  =  T;i  +Tjj+Tkk 

Where  the  above  subscripts  i,  j,  and  k  denote  the  corresponding  components  of  the  specified 
vector.  With  the  above  relations  for  the  components  of  the  thickness  vector,  the  perimeter’s 
contributions  to  presented  area  will  be  the  absolute  sum  of  the  k  components  of  the  cross 
products  of  L  and  T  and  W  and  T. 

4.2.2  Pitch.  Yaw,  and  Roll. 

A  three-dimensional  description  of  the  fragment’s  orientation  at  impact  has  been  included  with 
the  results.  Although  this  description  is  not  needed  for  the  penetration  equations,  it  may  expand 
the  utility  of  this  data  to  other  modeling  efforts  under  this  program. 

A  pitch,  yaw ,  and  roll  orientation  description  was  chosen.  This  orientation  was  chosen  as  it  is  a 
standard  input  parameter  for  other  empirical  models  which  may  be  of  potential  use  to  this 
program  [3].  Figure  6  illustrates  a  three-dimensional  fragment  in  space  prior  to  experiencing  a 
rotation.  The  rectangular  fragment  may  be  defined  in  terms  of  single  length,  width,  and  thickness 
vectors  (L,  W,  and  T).  To  ease  description  and  calculation  of  the  fragment’s  impact  orientation, 
unit  vectors  (Lu  =  L/Length,  Wu  =  W/Width,  and  Tu  =  T/Thickness)  will  be  used  to  describe 
the  transformed  axes  of  the  length,  width,  and  thickness  vectors. 
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FIGURE  6.  THREE-DIMENSIONAL  FRAGMENT  AND  ORIENTATION  UNIT  VECTORS 


The  fragment  experiences  three  separate  and  sequential  rotations  about  each  axis  as  illustrated  in 
figure  7.  Upon  completion  of  these  rotations,  a  description  of  the  fragment  orientation  may  be 
defined  by  the  angular  displacement  of  each  of  these  axis  rotations  (i.e.,  pitch,  yaw,  and  roll). 


FIGURE  7.  PITCH,  YAW,  AND  ROLL 

The  first  rotation  experienced  is  a  pitch  rotation  about  the  width  unit  vector  axis.  As  shown  in 
figure  7,  the  Wu  axis  does  not  move  during  this  rotation,  but  the  Lu  and  Tu  axes  are  translated 
an  equivalent  radial  pitch  distance.  It  should  be  noted  that  all  the  rotations  displayed  in  figure  7 
are  positive  rotations  by  the  right-hand  rule  method.  For  the  defined  pitch  axis,  a  downward 
pitch  angle  is  positive. 
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The  second  rotation  is  a  yaw  rotation  about  the  translated  Tu  axis.  The  Lu  and  Wu  axes  are 
rotated  an  equivalent  radial  yaw  distance  while  the  Tu  axis  does  not  move.  The  final  rotation  is 
a  roll  rotation  about  the  twice  translated  Lu  axis.  The  Lu  axis  does  not  move  during  this  rotation 
while  the  Tu  and  Wu  axes  are  rotated  an  equivalent  radial  roll  distance. 


It  is  reemphasized  that  the  three  rotations  described  above  must  be  conducted  in  the  specified 
sequence  (first  pitch,  then  yaw,  finally  roll).  Any  other  combination  of  rotations  may  result  in  a 
different  final  orientation  of  the  axes.  The  final  orientation  of  the  Wu,  Tu,  and  Lu  axes  to  their 
original  orthogonal  i,  j,  and  k  orientation  may  be  described  by  the  following  matrix  [3]: 


Note:  S  and  C  denote  sin  and  cos.  Pitch,  yaw,  and  roll  are  denoted  by  p,  y,  and  r. 


Lu 

Wu 

= 

Tu 

Sy 

CyCr 

-CySr 


-SpCy 

CpSr+SpSyCr 

CpCr-SpSySr 


CpCy 

SpSr-CpSyCr 

SpCr+CpSySr 


i 

j 

k 


Or  by  expansion  of  this  matrix: 


Lu  =  sin  (yaw)  i  -sin  (pitch)  cos  (yaw)  j  +  cos  (pitch)  cos  (yaw)  k 

Wu  =  cos  (yaw)  cos  (roll)  i  +  [cos  (pitch)  sin  (roll)  +  sin  (pitch)  sin  (yaw)  cos  (roll)]  j  + 
[sin  (pitch)  sin  (roll)  -  cos  (pitch)  sin  (yaw)  cos  (roll)]  k 

Tu  =  -cos  (yaw)  sin  (roll)  i  +  [cos  (pitch)  cos  (roll)  -  sin  (pitch)  s  (yaw)  sin  (roll)]  j  + 
[sin  (pitch)  cos  (roll)  +  cos  (pitch)  sin  (yaw)  sin  (roll)]  k 

The  above  relations  provide  ample  definition  of  the  orientation  of  a  rectangular  fragment  based 
on  its  length,  width,  and  thickness  vectors.  This  orientation  is  defined  through  application  of 

three  angular  translations  as  denoted  by  pitch,  yaw,  and  roll. 

5.  RESULTS. 

5.1  DATA  PRESENTATION. 

The  data  and  results  from  each  test  series  is  presented  in  appendix  B  in  both  tabular  and  graphic 
format.  These  include  a  tabulation  of  the  test  data  and  predictions  along  with  graphical 
comparisons  of  velocity  predictions,  prediction  accuracy,  and  median  comparison. 

5.1.1  Data  Tabulation. 

The  data  for  each  test  series  is  presented  in  tabulated  form.  Column  #  1  of  each  data  table  lists 
the  reference  shot  number.  Columns  #  2  and  3  lists  each  fragment’s  weight  and  dimensions. 
Columns  #  4  and  5  lists  the  fragment’s  impact  orientation  and  presented  area.  Column  #  6  and  8 
presents  the  fragment’s  impact  and  residual  velocity  as  measured  from  the  high-speed  film 
Columns  #  7  and  9  present  the  fragment’s  residual  velocity  and  ballistic  limit  velocity  as 
predicted  by  the  penetration  equations.  These  predictions  were  based  upon  the  target  and  the 
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fragment’s  characteristics  and  impact  conditions.  Finally,  column  #10  presents  the  residual 
velocity  prediction  accuracy  of  the  penetration  equations. 

5.1.2  V elocitv  Comparisons . 

For  each  test  series,  a  direct  comparison  of  the  predicted  to  the  actual  residual  velocities  is 
presented  graphically  in  columnar  format.  The  graphs  also  include  impact  velocities  for  each 
comparison  as  an  aid  in  visually  interpreting  the  magnitude  of  prediction  error. 

5.1.3  Prediction  Accuracy. 

A  second  graph  presenting  the  prediction  accuracy  of  the  penetration  equations  is.  included  for 
each  test  series.  The  relation  utilized  is  as  follows: 

Prediction  Accuracy  =  1  -  (Vr  -  Vr,UDM)/Vl 

The  relation  measures  the  difference  between  the  predicted  and  actual  residual  velocity  with 
respect  to  the  fragment’s  impact  velocity.  In  this  manner,  the  magnitude  of  the  difference 
between  the  predicted  and  actual  residual  velocities  may  be  kept  in  perspective  to  the  original 
impact  velocity.  An  accuracy  value  close  to  1  represents  an  accurate  prediction.  Accuracy 
values  greater  than  1  represent  residual  velocity  predictions  which  are  greater  than  the  actual 
residual  velocity.  Also,  accuracy  values  less  than  1  represent  residual  velocity  predictions  which 
are  less  than  the  actual  residual  velocity. 

The  range  of  prediction  accuracy  varies  from  zero  to  2.  A  prediction  accuracy  of  zero  would 
represent  a  prediction  of  nonpenetration  while  in  actuality  the  residual  velocity  equals  the  impact 
velocity.  An  accuracy  of  2  would  represent  the  opposite  extreme,  with  a  predicted  residual 
velocity  equal  to  the  impact  velocity,  while  the  actual  residual  velocity  equates  to  zero. 

5.1.4  Median  Comparison. 

Residual  velocity  vs.  impact  velocity  prediction  plots  are  presented  for  each  test  series.  The 
actual  residual  velocity  data  points  have  also  been  plotted  for  comparison  with  the  prediction 
plots.  The  fragmentation  characteristics  utilized  in  the  predictions  were  the  average  fragment 
weight  and  dimensions  for  each  test  series. 

Each  graph  includes  plots  for  the  possible  orientation  extremes  of  the  fragment.  These  include 
plots  for  both  a  completely  flat  and  an  edge-on  impact.  An  additional  plot  is  included  which 
represents  fragment  impacts  occurring  at  the  median  presented  area.  The  value  for  this  median 
area  was  calculated  from  penetrating  fragments  only. 

The  median  presented  area  was  chosen  as  this  would  allow  a  means  for  direct  comparison  with 
the  data  points.  As  such,  50%  of  the  data  points  that  display  residual  velocity  (i.e.,  penetrating 
fragments)  would  be  expected  to  be  above  the  median  prediction  line.  The  other  50%  would  be 
expected  to  be  below. 
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5^2.  LIMITATIONS  OF  THE  PENETRATION  EQUATIONS 

There  are  many  factors  which  may  account  for  the  limitations  in  prediction  accuracies  of  the 
penetration  equations.  A  discussion  follows  which  outlines  factors  which  were  of  the  most 
significance  to  the  accuracy  of  the  equations.  Additional  factors  covered  include  those  which 

may  intuitively  appear  to  be  of  significance  but  through  further  analysis  were  eliminated  from 
consideration. 

5.2.1  Ballistic  Perforation  Mode. 

The  ballistic  limit  velocity  (V50)  is  an  important  input  variable  for  the  residual  velocity  equation. 
An  error  for  this  value  may  result  in  significant  error  in  the  predicted  value  for  residual  velocity. 
The  value  for  V50  is  calculated  through  the  use  of  the  empirically  derived  ballistic  limit  equation. 

The  ballistic  limit  equation  assumes  that  the  fragment’s  presented  area  is  the  only  influential 
fragment  shape  factor  affecting  perforation.  It  can  be  easily  understood  that  this  may  be  an 
madequate  assumption.  As  an  example,  consider  the  two  penetrators  of  figure  8.  They  both  are 
impacting  a  target  at  the  same  presented  area.  They  are  of  the  same  mass  and  material,  but  they 
differ  in  impact  shape.  One  is  pointed  while  the  other  is  blunt.  The  ballistic  limit  equation 
would  predict  an  equivalent  ballistic  limit  for  each  impact.  However,  we  would  intuitively 
expect  these  values  to  be  different. 


FIGURE  8.  EXAMPLE  POINTED  AND  BLUNT  PENETRATORS 

The  data  points  which  were  used  in  Project  Thor  [4]  to  develop  the  empirical  constants  for  the 
ballistic  limit  equation  were  accumulated  through  the  use  of  compact  fragment  simulating 
projectiles  (FSP)  [5].  The  geometry  of  these  FSP  were  cubes  and  cylinders  with  length  to 
diameter  ratios  close  to  one  [6],  A  relatively  blunt  impact  was  thus  assured  regardless  of  the 
orientation  of  the  FSP  to  the  target.  As  such,  it  may  be  assumed  that  plugging  was  the 
predominant  perforation  mode  observed  during  the  development  of  the  ballistic  limit  equation. 
However,  very  few  of  the  test  shots  outlined  in  this  report  resulted  in  a  completely  sheared  plug 
which  separated  from  the  target.  Figure  9  outlines  the  two  predominant  perforation  modes 
observed  during  these  tests.  These  were  petaling  and  plugging. 
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Petaling 
(Single  Sided) 


FIGURE  9.  OBSERVED  PERFORATION  MODES 


5.2.2  Petaling. 

In  petaling  perforation,  a  sharp-edged  fragment  creates  an  intense  shear  stress  over  a  small  region 
of  the  target.  This  results  in  an  initial  breech  of  the  target  material  near  the  impact  area.  This 
material  remains  attached  to  the  plate  but  accelerates  outward  incurring  an  initial  velocity  and 
thus  kinetic  energy  from  the  penetrator  impact.  As  the  material  is  pushed  outward,  it  continues 
its  acceleration  (and  kinetic  energy  rise)  but  swings  away  from  the  path  of  the  penetrator, 
allowing  relatively  easy  passage  of  the  fragment.  However,  postimpact  collision  with  the 
material  can  occur  and  may  result  in  a  deflection  of  the  fragment  from  its  original  trajectory. 

As  the  plate  material  remains  partially  attached  to  the  plate,  additional  shearing  of  the  target  may 
occur  as  the  material  decelerates  and  dissipates  its  kinetic  energy.  Dissipation  of  the  kinetic 
energy  also  occurs  through  deformation,  resulting  in  a  curved  surface  which  resembles  a  flower 
petal  (i.e.,  petaling).  For  thin  target  plate  material,  the  damage  to  the  target  from  shearing  may  be 
similar  to  that  of  tearing  paper.  The  resultant  perforation  hole  may  be  much  larger  than  the 
original  impact  presented  area.  It  is  important  to  stress  that  the  majority  of  these  deformations 
occur  after  the  passage  of  the  fragment.  As  such,  they  have  little  influence  on  the  fragment. 

5.2.3  Plugging. 

Plugging  occurs  primarily  for  blunt  edged  penetrator  impacts.  For  thin,  ductile  target  materials, 
the  target  may  bulge  or  dish  out  around  the  perimeter  of  the  penetrator’s  impact  surface. 
Eventually,  shearing  stresses  become  so  great  along  this  perimeter  that  failure  of  the  plate 
material  occurs  there.  This  newly  formed  “plug”  of  plate  material  remains  in  the  path  of  the 
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projectile.  It  is  thus  accelerated  forward,  resulting  in  further  deceleration  of  the  fragment.  These 

accelerations  are  completed  when  the  fragment  and  plug  eventually  reach  an  equivalent  residual 
velocity. 

Figures  10  and  11  provide  excellent  examples  of  petaling.  Both  figures  are  target  rear  views  of 
impact  damage.  Figure  10  resulted  from  the  impact  of  a  fan  blade  fragment.  Although  the  fan 
blade  fragment  was  relatively  large,  the  resulting  hole  and  “petal”  is  much  larger  than  the 
presented  area  and  the  actual  size  of  the  fragment.  Figure  11  resulted  from  a  turbine  blade 
impact  and  illustrates  that  more  than  one  petal  may  form. 


FIGURE  10.  PETALING  (SINGLE  SIDED)  FIGURE  11.  PETALING  (MULTISIDED) 

Figure  12  resulted  from  impact  of  a  simulated  blade  fragment  from  test  series  2.  A  plug  of  metal 
has  been  completely  sheared  from  the  target.  The  area  which  was  sheared  is  approximately  equal 
to  that  of  the  presented  area  of  the  fragment.  As  this  fragment  impacted  relatively  flat,  the  entire 
outline  of  the  hole  is  approximately  equal  to  that  of  the  fragment’s  side.  Additionally,  note  the 
dishing  out  of  the  target  surface  near  the  edges  of  the  hole. 

Finally,  figure  1 3  has  been  included  to  emphasize  that  a  combination  of  the  discussed  perforation 
modes  may  occur.  The  figure  provides  illustration  of  plugging  and  petaling.  The  perforation 
first  started  as  plugging.  Note  the  dishing  out  of  the  plate  material,  along  with  the  nice  outline  of 
the  fragment  near  the  top  edges  of  the  hole.  However,  a  portion  of  this  plug  remained  affixed  to 
the  target  after  penetration.  As  a  result,  a  petal  has  formed.  The  imprint  of  the  entire  fragment’s 
side  can  be  clearly  seen  on  the  edge  of  this  “petal.”  Additionally,  the  resulting  damage  hole  is 

noted  as  being  much  larger  than  that  of  the  fragment’s  imprint. 
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FIGURE  12.  PLUGGING 


FIGURE  13.  PLUGGING  AND  PETALING 


5.2.4  Ballistic  Perforation  Mode  and  Fragment  Orientation. 

At  the  representative  velocities  of  an  uncontained  engine  blade  fragment,  the  main  factors  which 
would  determine  the  type  of  ballistic  perforation  mode  encountered  are  the  geometry  and 
orientation  of  the  fragment.  Figure  14  illustrates  three  fragments  at  various  orientation  angles. 
Note  that  orientation  angle  is  defined  as  the  angle  between  the  perpendicular  of  the  fragment’s 
trajectory  to  the  plane  of  the  fragment. 


0  =  0* 


Fragment 

Trajectory 


In  general,  plugging  perforation  occurs  mainly  for  blunt  projectile  impacts.  For  the  sharp-edged 
engine  blade  fragments,  this  would  occur  only  for  fragment  orientation  angles  close  to  zero 
degrees.  At  orientation  angles  close  to  90  degrees,  we  would  expect  multisided  petaling  as  the 
target  shears  along  both  sides  of  the  sharp-edged  fragment.  At  orientation  angles  between  0°  and 
90°,  we  would  expect  single-sided  petaling.  There  are  some  exceptions,  such  as  impacts 
occurring  along  the  flat  edges  of  the  simulated  fragment. 


Through  analysis  of  the  test  data,  it  was  determined  that  the  mean  fragment  presented  area  of  all 
tests  was  approximately  50%  of  a  fragment’s  maximum  side  area.  This  equates  to  a  mean 
orientation  angle  of  60  degrees.  Thus  we  would  expect  that  single-sided  petaling  would  be  the 
predominant  perforation  mode  encountered  during  testing.  This  is  confirmed  by  examination  of 
the  target  photographs  of  appendix  C. 

—  2-5  Comparison  of  Perforation  Mode  to  Penetration  Efficiency. 

The  penetration  efficiency,  or  ease  of  penetration  for  each  of  the  two  predominant  perforation 
modes,  may  vary  with  respect  to  each  other  [5],  For  a  specified  fragment,  the  variance  in  these 
efficiencies  is  mainly  dependent  on  the  thickness  and  material  properties  of  the  target  material. 
The  following  discussion  assumes  a  relatively  thin  and  ductile  target  material,  such  as  aircraft 
aluminum. 

Consider  again  the  fragments  of  figure  9.  Each  are  of  the  same  mass  and  material.  Additionally, 
they  each  project  the  same  presented  area  to  the  target.  However;  they  differ  in  shape,  and 
orientation  with  respect  to  the  target.  They  may  also  differ  in  impact  velocity  such  that 
penetration  of  their  respected  target  is  assured. 

The  fragment  which  resulted  in  petaling  passed  through  the  target  with  much  more  ease  than  the 
plugging  fragment.  This  is  due  to  the  shear  concentrations  which  formed  at  the  interface  of  the 
petal  with  the  target.  Also,  as  the  petal  continued  formation,  the  plate  material  moved  away  from 
the  direction  of  the  fragment’s  path.  This  effect  is  especially  pronounced  for  single  sided 
petaling.  Additionally,  the  energy  intensive  dishing  deformation  observed  during  plugging  is 
absent  in  petaling.  In  general,  it  may  be  stated  that  for  aircraft  aluminum  targets,  petaling  is  a 
more  efficient  (less  energy  absorbing)  perforation  mode  than  plugging. 


The  V50  values  for  petaling  and  plugging  may  be  ranked  with  the  least  efficient  having  the 
highest  ballistic  limit  value.  As  such,  for  a  specified  fragment  at  a  specified  impact  orientation,  it 
may  be  assumed  that 


V50, Plugging  >  V50,Petaling 

5.2.6  Perforation  Mode  and  Prediction  Accuracy. 


Recall  that  the  ballistic  limit  equation  provides  an  empirical  estimate  for  V50.  The  penetration 
data  which  was  utilized  to  develop  this  equation  was  mainly  based  on  the  plugging  type  of 
perforation  mode.  Therefore,  it  may  be  assumed  that  for  the  input  conditions  utilized  by  this 
equation,  it  may  provide  a  fair  estimate  of  V50  should  plugging  occur.  Now  consider  the  residual 
velocity  equation 


J 


V2_y2 


50. 


PA  t 

1.0+ - tL- 


Wcos© 


The  Residual  Velocity  Equation 


18 


If  plugging  is  predominant,  the  equation  should  provide  a  good  estimate  of  residual  velocity  as 
the  input  value  for  V50  should  be  reasonably  accurate.  But  what  if  petaling  is  the  predominant 
perforation  mode?  As  previously  discussed,  the  ballistic  limit  velocity  for  petaling  would  be 
expected  to  be  lower  than  for  plugging  (and  as  such,  lower  than  that  predicted  by  the  ballistic 
limit  equation).  Thus  by  the  residual  velocity  equation,  should  petaling  be  the  observed 
perforation  mode,  we  would  expect  the  actual  residual  velocity  to  be  greater  than  predicted 

VsO,Predicted  £  V 50, Plugging ,  then  V 50,Petaling<  Vso^redicted 

therefore  (V2  -  V25oJpetaiing)  >(V2  -  V250>Predicted),  or  Vr,  Petaling  ^Predicted 

This  can  be  confirmed  by  comparing  the  prediction  accuracy  of  each  test  shot  with  its  respective 
photograph  in  appendix  C.  We  note  in  these  photographs  a  predominance  of  petaling  damage  for 
shots  with  actual  Vr  »  predicted  Vr. 

Another  causative  factor  may  be  that  plugging  is  a  major  assumption  which  was  utilized  in  the 
derivation  of  the  residual  velocity  equation.  It  is  assumed  that  at  impact,  a  plug  of  the  target 
material  is  sheared  from  the  target.  The  plug  is  also  assumed  to  be  equal  in  shape  and  area  to  that 
of  the  presented  area  of  the  fragment.  After  formation,  the  plug  is  assumed  to  depart  the  target  at 
the  same  velocity  as  the  residual  velocity  of  the  fragment. 

Another  assumption  utilized  in  the  derivation  is  that  after  penetration  the  fragment  continues  to 
travel  along  its  original  trajectory.  This  assumption  is  important  for  portions  of  the  derivation 
dealing  with  the  conservation  of  momentum.  In  reality,  deflection  of  the  fragment  from  its 
original  trajectory  is  possible.  As  these  deflections  are  not  accounted  for,  prediction  inaccuracies 
may  be  expected. 

5.2.7  Perimeter  of  Shear. 

The  empirical  constants  of  the  ballistic  limit  equation  were  developed  under  the  assumption  that 
presented  area  was  the  main  fragment  shape  factor  affecting  perforation.  However,  assuming  that 
plugging  perforation  is  the  predominant  failure  mode,  the  actual  shearing  and  failure  of  the  target 
occurs  at  the  perimeter  of  the  perforation  plug.  As  such,  it  is  suggested  that  the  ballistic  limit 
may  also  have  some  dependence  on  the  perimeter  of  the  presented  area. 

The  fragments  utilized  in  the  development  of  the  empirical  constants  were  cubic  and  cylindrical 
in  shape.  The  presented  area  and  perimeter  of  the  presented  area  for  these  shapes  do  not  change 
drastically  with  the  orientation  of  these  fragments.  However,  engine  blade  fragments  are  in 
general  relatively  thin  rectangles.  As  such,  both  the  presented  area  and  the  perimeter  of  the 
presented  area  may  vary  greatly  with  orientation  of  these  fragments. 

Consider  the  three  fragments  and  presented  areas  of  figure  15.  All  three  presented  areas  are 
equivalent.  Also,  all  three  fragments  are  of  the  same  mass  and  material.  They  only  differ  in  the 
shape  of  presented  area  as  rectangles  of  decreasing  width  to  length  ratios. 
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FIGURE  15.  PERIMETER  COMPARISON 

For  a  given  target,  the  ballistic  limit  equation  would  predict  the  same  V50  value  for  all  three 
fragments,  as  their  weights  and  presented  areas  are  equivalent.  But  for  each  fragment,  consider 
the  actual  length  of  shear  damage  which  would  occur  by  plugging  perforation.  We  would  expect 
the  extent  of  this  damage  to  be  proportional  to  the  energy  expended  by  perforation.  If  this 
shearing  occurs  at  the  perimeter  of  the  presented  area  (as  defined  by  plugging),  then  the  expended 
energy  should  also  be  proportional  to  these  perimeters. 

Figure  16  compares  the  perimeter  of  a  rectangle  to  the  perimeter  of  a  square  of  equivalent 
presented  area.  The  plot  illustrates  how  the  perimeter  of  the  rectangle  increases  with  respect  to 
the  square  as  the  rectangles  width  to  length  ratio  decreases.  Note  that  for  width  to  length  ratios 
above  0.6,  the  perimeter  of  the  square  and  rectangle  are  approximately  equal.  But  for  width  to 
length  ratio  values  below  0.4,  the  rectangle’s  perimeter  becomes  significantly  greater  than  the 
square  s  perimeter.  Figure  17  is  included  as  supportive  evidence  that  the  perimeter  ratio  does 
have  an  effect  on  the  accuracy  of  the  predictions. 

As  displayed  in  figure  17,  individual  test  prediction  accuracies  have  been  plotted  against  the  ratio 
of  each  fragment  s  presented  area  to  its  maximum  side  area.  This  ratio  is  approximately  equal  to 
the  presented  area’s  rectangular  width  to  length  ratio  for  fragments  with  side  width  to  length 
ratios  of  greater  than  0.6  (all  except  the  turbine  blades).  The  side  area  of  fragments  in  this  range 
approximates  that  of  a  square. 
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FIGURE  16.  PERIMETER  RATIO  (RECTANGLE  AND  SQUARE  OF  EQUIVALENT 

PRESENTED  AREA) 
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FIGURE  17.  PREDICTION  ACCURACY  VERSUS  AREA  RATIO 
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Note  in  figure  17  the  predominance  of  underpredictions  (with  accuracies  less  than  one)  at  ratios 
greater  than  0.6.  In  this  range,  the  perimeter  of  the  roughly  square  presented  area  is  in  agreement 
with  the  penetration  equation’s  coefficients.  The  majority  of  underpredictions  noted  in  this  range 
are  as  a  result  of  petaling.  However,  note  that  as  the  ratio  falls  below  0.4,  prediction  accuracies 
approach  one.  This  is  not  that  the  penetration  equations  become  more  accurate  in  this  range,  but 
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rather  that  the  effects  of  petaling  are  being  masked  by  increasing  overprediction  errors  caused  by 
the  increasing  perimeter  of  the  presented  area.  Eventually  as  the  ratio  falls  even  lower,  the 
predictions  all  become  overpredictions  as  the  perimeter  effect  becomes  pronounced  at  near  edge- 
on  impacts. 

5,2,8  Presented  Area  Accuracy. 

Presented  area  of  the  fragment  to  the  target  is  an  important  dependent  variable  in  both  the 
residual  velocity  and  ballistic  limit  equations.  The  accuracy  of  this  input  variable  is  essential  when 
impact  conditions  are  close  to  the  ballistic  limit.  Presented  area  accuracy  is  also  critical  when  the 
ratio  of  the  fragment’s  maximum  to  minimum  possible  presented  area  is  large  (i.e.,  a  large,  thin, 
flat  fragment).  This  may  be  seen  in  figures  18  and  19  for  the  simulated  fragment  and  compressor 
blade  respectively. 

The  presented  area  ratio  for  the  simulated  fragment  is  relatively  small  (approximately  5.6).  In 
figure  18,  the  bandwidth  (or  distance)  between  edge-on  and  flat  impact  predictions  is  also 
relatively  small.  However,  for  the  compressor  blade  the  presented  area  ratio  is  much  larger 
(approximately  28.6)  than  for  the  simulated  fragment.  As  such,  in  figure  19  the  bandwidth 
between  edge-on  and  flat  impact  predictions  is  also  much  larger.  Due  to  this  large  bandwidth, 
ballistic  predictions  for  the  compressor  blade  are  much  more  sensitive  to  presented  area  input 
accuracy  than  for  the  simulated  fragment. 


UDM  Edge  Impact  . ""  UDM  Flat  Impact  '  UDM  Median  Penetrators 


FIGURE  18.  TEST  SERIES  2,  RESIDUAL  VELOCITY  CURVES  ORIENTATION 

COMPARISON 
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FIGURE  19.  TEST  SERIES  5  AND  7,  RESIDUAL  VELOCITY  CURVES  ORIENTATION 

COMPARISON 

Figure  20  (from  test  series  2)  illustrates  the  importance  of  presented  area  input  accuracy  for 
impact  conditions  which  are  close  to  the  ballistic  limit.  In  this  figure,  the  presented  area  is  varied 
by  a  relatively  small  amount  of  ±  10%.  In  shots  14,  31,  and  32,  the  predictions  are  only  minorly 
affected  by  this  variation. 
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500 

450 

CO  400 
i£  350 

>  300 

O  250 

O 

o  200 

>  150 
100 

50 

0 


FIGURE  20.  EFFECTS  OF  PRESENTED  AREA  INACCURACIES 
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The  variations  resulted  in  dramatic  prediction  differences  for  shots  33  and  34.  An  increase  in 
presented  area  resulted  in  predictions  of  nonpenetration  for  both  shots  (or  residual  velocity  =  0, 
as  denoted  by  the  absence  of  the  respective  bar  plot).  As  such,  we  would  expect  these  impacts  to 
be  relatively  close  to  their  ballistic  limit.  This  is  confirmed  by  comparing  the  impact  velocities 
of  these  shots  to  their  predicted  ballistic  limits.  For  shot  33,  the  impact  velocity  was  416  fps  with 
a  predicted  ballistic  limit  of  384  fps.  The  impact  velocity  of  shot  34  was  375  fps  with  a  predicted 
ballistic  limit  of  350  fps. 

As  discussed,  presented  area  must  be  accurately  inputted  to  the  penetration  equations  to  assure 
accurate  predictions.  Unfortunately,  the  presented  area  must  be  calculated  from  visual  high¬ 
speed  film  inputs  which  are  measured  during  the  dynamic  event  of  impact.  There  are  numerous 
factors  which  affect  the  accuracy  of  these  calculations.  This  includes  poor  visibility  of  the 
fragment  at  impact.  Additionally,  the  fragment’s  orientation  (and  thus  presented  area)  is 
measured  one  frame  prior  to  impact  to  negate  the  possibility  of  measuring  orientation 
inaccuracies  from  impact  induced  rotation.  However,  if  the  fragment  is  tumbling,  some  rotation 
of  the  fragment  will  occur  between  film  frames,  resulting  in  further  presented  area  inaccuracies. 
Finally,  the  fragments  were  assumed  to  be  rectangular  in  shape  to  ease  calculation  of  the 
presented  area.  In  reality,  the  fragments  are  somewhat  differently  shaped,  owing  to  their  airfoil 
characteristics.  The  compressor  blades  are  a  notable  example,  with  an  aerodynamic  twist  of 
approximately  30  degrees  along  their  length. 

5.2.9  Target  Material  and  Structural  Characteristics. 

It  is  assumed  that  new  metal  targets  were  utilized  in  the  development  of  the  empirical  constants 
utilized  by  the  ballistic  limit  equation.  Additionally,  the  way  these  targets  were  supported  during 
these  tests  is  unknown.  Both  of  these  may  be  important  factors  in  the  development  of  these 
constants. 

New  2024  T3  aluminum  targets  were  utilized  during  test  series  2  through  7.  For  most  of  these 
tests,  the  fragment  impact  points  were  distant  from  the  target  support  structure.  For  these  test 
series,  predictions  of  the  residual  velocity  were  noted  to  be  relatively  close  to  that  of  the  actual 
measured  residual  velocity.  The  notable  exception  to  this  were  the  results  from  test  series 
number  8  which  impacted  the  actual  aluminum  aircraft  structure. 

Test  series  8  utilized  P-3  Orion  aircraft  engine  cowlings  as  impact  targets  (see  figures  21  and  22). 
These  cowlings  are  shaped  and  supported  in  numerous  spots  by  aircraft  stringers.  Additionally, 
during  their  normal  lifetime,  the  cowlings  are  expected  to  have  been  exposed  to  many  cyclic 
extremes  in  atmospheric  temperature  and  pressure. 

Figure  23  illustrates  the  velocity  comparisons  for  test  series  8.  Of  interest  is  the  stringer  impacts 
which  occurred  for  approximately  half  of  the  test  impacts.  It  should  be  noted  that  the  intended 
impact  points  were  midway  between  the  stringers  for  each  shot.  Fragment  tumbling  resulted  in 
apparently  random  trajectory  deflections  causing  the  stringer  impacts. 
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FIGURE  21.  COWLING  EXTERIOR 


FIGURE  22.  COWLING  INTERIOR 
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FIGURE  23.  TEST  SERIES  8  PREDICTION  COMPARISON 

Due  to  the  randomness  of  the  trajectory  deflections,  it  may  be  possible  to  predict  the  occurrence 
of  these  impacts.  This  could  be  achieved  by  randomly  plotting  (on  a  computer)  the  previously 
mentioned  mean  presented  area  (50%  of  the  maximum  side  area)  within  the  area  bounding  the 
stringers.  Instances  where  the  presented  area  overlaps  a  stinger  would  be  counted  and  compared 
with  the  total  number  of  plots  to  arrive  at  a  probability  of  stringer  impact. 
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The  residual  velocities  for  tests  which  resulted  in  stringer  impacts  are,  as  intuitively  expected, 
much  lower  than  predicted.  As  such,  it  may  be  concluded  that  stringer  placement  density  does 
provide  some  mitigation  of  uncontained  engine  failure. 

Figure  24  illustrates  velocity  predictions  for  cowling  tests  which  resulted  in  non-stringer  impacts. 
Of  interest  is  the  large  underprediction  of  residual  velocity  for  shots  60  and  67.  As  seen  in  the 
impact  photographs  of  figures  25  and  26,  these  shots  did  result  in  actual  plugging  failure.  We 
would  expect  the  predictions  for  these  tests  to  be  relatively  accurate. 


o  V  impact  0  V  res  a  V  res  UDM 


Shot  Number 

FIGURE  24.  TEST  SERIES  8  PREDICTIONS  (SKIN  IMPACTS  ONLY) 


(a)  (b) 

FIGURES  26.  PLUGGING  FAILURE,  SHOT  67  FRONT  (a)  AND  BACK  (b) 

Explanation  of  these  results  are  difficult.  One  causative  factor  may  be  due  to  the  unavailability 
of  penetration  equation  coefficients  for  the  material  that  the  cowlings  were  composed  of  (7075 
T6  Aluminum).  Coefficients  for  2024  T3  aluminum  were  used  instead.  Although  7075  is  of 
increased  strength  than  2024,  it  is  also  more  brittle.  As  such,  it  may  be  more  susceptible  to 
impact  failures. 

Another  explanation  may  attribute  these  results  to  changes  in  the  material  characteristics  of  the 
cowlings  which  were  incurred  during  their  normal  lifetime  of  service.  Environmental  extremes 
which  these  cowlings  experienced  may  have  resulted  in  a  tempering  of  the  aluminum,  making  the 
metal  more  susceptible  to  impact. 

Another  explanation  may  be  that  in  addition  to  temperature  extremes,  these  cowlings  were 
exposed  to  high  degrees  of  humidity  in  a  salt  laden  atmosphere.  The  P-3  is  utilized  as  a  maritime 
patrol  aircraft  and  normally  flies  missions  at  very  low  altitudes  over  the  ocean  for  extended 
durations.  Corrosion  of  the  metal  from  numerous  exposures  to  the  sea  air  may  have  weakened 
the  ballistic  resistance  of  these  targets. 

Finally,  all  of  the  impacts  occurred  very  close  to  or  on  top  of  an  aircraft  stringer.  This  is  due  to 
the  high  placement  density  of  stringers  on  this  structure.  It  is  clear  that  those  fragments  which 
impacted  a  stringer  would  probably  have  a  lower  residual  velocity  than  predicted.  However,  it  is 
unclear  as  to  how  the  stringers  may  contribute  to  impacts  which  occur  very  near,  but  not  on  top 
of,  their  placement. 

5.2.10  Rotational  Versus  Linear  Kinetic  Energy. 

As  previously  discussed,  many  of  the  fragments  tumbled  severely  due  to  their  aerodynamic 
shape.  Some  concern  may  arise  that  the  influence  of  the  fragment’s  rotational  kinetic  energy  may 


be  of  significance  in  the  ballistic  predictions.  The  most  significant  tumble  occurred  for  tests 
which  utilized  the  compressor  blades.  The  aerodynamic  twist  of  these  blades  induced  a  high  rate 
of  rotation  along  their  longitudinal  axis. 


The  kinetic  energy  (rotational  and  linear)  may  be  calculated  by  the  following  relation: 

KErotational  =  I© 

KE, inear  =  1/2  mV2 


I  =  Mass  moment  of  inertia  about  the  rotational  axis 
©  =  Rotational  velocity  (radians/s) 
m  =  Fragment’s  mass 
V  =  Impact  velocity 

Assuming  a  thin  rectangular  plate  as  the  fragment,  the  mass  moment  of  inertia  about  the 
fragment’s  longitudinal  axis  may  be  calculated  as  I  =  1/12  mW2,  where  W  is  the  fragment’s 
width  (for  the  compressor  blades,  W  =  2.75  inches). 

Taking  care  to  utilize  consistent  units,  a  comparison  of  the  rotational  to  kinetic  energy  for  the 
compressor  blade  fragments  is  presented  in  table  1  and  displayed  graphically  in  figure  27.  These 
presentations  clearly  illustrate  that  although  the  fragments  are  tumbling  at  a  very  high  rotational 
rate,  their  rotational  kinetic  energy  is  negligible  in  comparison  to  their  linear  (or  translational) 
kinetic  energy.  As  such,  the  influence  of  the  rotational  kinetic  energy  is  also  negligible  to  the 
dynamics  of  impact. 


TABLE  1.  TEST  SERIES  5  AND  7,  KINETIC  ENERGY  TABULAR  COMPARISON 


Shot 

Number 

Fragment 

Mass 

(grams) 

Rotational 
Velocity  Along 
Longitudinal  Axis 
(revs/s) 

Linear  Impact 
Velocity 
(ft/s) 

Rotational 

Impact 

Kinetic  Energy 
(fUbs) 

Linear  Kinetic 
Energy 
(ft.lbs) 

41 

71.9 

165 

539 

23 

716 

42 

70.9 

99 

508 

8 

627 

43 

71.5 

62 

366 

3 

328 

45 

73.1 

124 

315 

13 

249 

46 

69.4 

99 

199 

8 

94 

47 

69.3 

52 

246 

2 

144 

48 

69.6 

26 

139 

1 

46 

49 

71.9 

41 

219 

1 

118 

50 

70.2 

50 

214 

2 

110 

51 

69.4 

0 

243 

0 

140 

52 

69.3 

71 

323 

4 

248 

53 

69.4 

71 

372 

4 

329 

54 

71.9 

0 

437 

0 

470 

55 

70.2 

99 

481 

8 

556 

56 

71.9 

83 

593 

6 

866 

57 

69.4 

83 

578 

6 

794 

58 

69.4 

0 

526 

0 

658 
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FIGURE  27.  TEST  SERIES  5  AND  7,  KINETIC  ENERGY  GRAPHICAL  COMPARISON 
5  3  DISCUSSION  OF  RESULTS. 

5,3,1  Prediction  Accuracy  (All  Test  Series). 

Figure  28  displays  the  distribution  histogram  based  on  the  prediction  accuracies  for  all  test  series. 
Unintentional  impacts  with  structural  components  (stringers)  which  occurred  during  test  series  8 
(cowling  tests)  are  not  included  within  the  histogram.  The  graph  displays  the  number  of 
predictions  which  fall  within  each  accuracy  range. 


_ _ 1 

I _ 

-- 

:  1 

_ 

i 

■ 

1  _ I 

n _ B 

1 _ _ _ 

i 

I 

3 _  _ _ E 

1 

Hi 

»-+ 

-H- 

U 

1,  ,[ 

n 

L-+ 

I 

i-H- 

1 — h 

— h 

r+ 

Shot  Number 


<0 

c 

o 


*o 

2 

CL 

«►- 

O 

%m 

O) 

.a 

E 

3 


Prediction  Accuracy 


FIGURE  28.  PREDICTION  ACCURACY  DISTRIBUTION 


29 


As  seen  in  the  graph,  the  majority  of  prediction  accuracies  were  less  than  one.  The  arithmetic 
mean  for  all  test  predictions  was  0.881  with  a  standard  deviation  of  0.36.  (Or  the  average 
prediction  fell  within  the  range  of  0.521  to  1.241).  It  may  be  concluded  that  for  the  majority  of 
the  test  shots,  the  penetration  equations  had  a  tendency  to  underpredict  the  residual  velocities  of 
the  fragments  (as  expected  due  to  petaling  induced  errors). 

5.3.2  Prediction  Accuracy  by  Fragment  Type. 

Figure  29  graphically  displays  the  mean  prediction  accuracy  for  each  type  of  fragment  tested.  The 
standard  deviation  for  each  of  these  accuracies  is  also  displayed. 
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FIGURE  29.  PREDICTION  ACCURACY  BY  FRAGMENT  TYPE 


Note  that  the  turbine  blades  provided  the  most  accurate  predictions  while  also  maintaining  the 
least  standard  deviation  (Mean  0.99/Standard  Deviation  0.14).  This  may  be  attributed  to  the 
highly  rectangular  shape  of  these  type  of  fragments.  Perimeter  shear  error  may  be  an  influential 
factor  for  these  type  of  fragments,  effectively  countering  target  petaling  induced  errors. 


The  least  accurate  predictions  were  those  calculated  for  the  engine  cowling  shots  with  a  large 
emphasis  of  underpredictions  of  residual  velocity  (Mean  0.65/Standard  Deviation  0.45).  It  should 
be  noted  that  those  tests  which  resulted  in  stringer  impacts  were  not  considered  in 
figure  29.  The  large  error  in  predictions  may  be  attributed  to  the  type  and  condition  of  the  target 
material.  The  only  ballistic  limit  empirical  constants  which  are  available  for  aluminum  are  those 
for  2024  T3  aluminum.  The  cowlings  were  constructed  of  7075  T6  aluminum.  As  such, 
2024  T3  aluminum  constants  were  utilized  as  a  substitute.  Also,  the  cowlings  had  been  aged 
through  a  normal  aircraft  life  cycle  of  environmental  extremes.  Additional  degradation  in  material 
condition  may  have  occurred  due  to  extended  open  air  storage  near  a  high  salt  content  dry  lake 
bed  (China  Lake). 
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Predictions  for  the  remainder  of  engine  fragments  (fan  and  compressor  blades)  were  similar  in 
accuracy  and  standard  deviation.  Mean  prediction  accuracies  for  these  fragments  were  fair 
(0.80-0.90)  with  the  predominance  of  underpredictions  being  attributed  to  petaling  induced  error. 
However,  standard  deviation  values  were  also  somewhat  large  for  these  fragments 
(0.26-0.39).  These  large  deviations  may  be  attributed  to  the  large,  but  thin  shape  of  the 
fragments.  Variations  in  fragment  orientation  for  each  impact  resulted  in  large  variations  in  the 
shape  of  the  presented  area.  As  such,  perimeter  shear  error  varied  greatly  between  each  shot, 
with  resultant  large  deviations  in  prediction  accuracy. 

5.3.3  Prediction  Accuracy  by  Prediction  Model. 

Large  deviations  in  prediction  accuracy  were  noted  for  predictions  of  fan  and  compressor  blade’s 
residual  velocities.  Several  alternative  prediction  models  for  the  ballistic  limit  were  examined  to 
determine  their  benefits  over  the  current  model.  These  models  were  evaluated  and  compared  in 
performance  with  the  JTCG/ME  ballistic  limit  equation  for  fragments. 

There  were  three  alternative  ballistic  limit  models  considered.  These  included  the  JTCG/ME 
ballistic  limit  equation  for  projectiles  [1].  This  equation  is  similar  to  the  fragment  equation,  but 
was  developed  to  better  account  for  petaling  effects.  Other  equations  considered  included  the 
“FAA  Energy  Equation”  [6].  This  equation  was  selected  for  comparison  due  to  its  inclusion  of 
perimeter  shear  as  a  dependent  variable.  Also,  a  modified  form  of  the  FAA  energy  equation  was 
evaluated  to  determine  prediction  enhancements  through  modification  of  the  equation’s  empirical 
constant. 


A  graphical  comparison  of  the  performance  of  these  models  is  displayed  in  figure  30.  The  chart 
outlines  the  mean  prediction  accuracy  of  each  model  as  a  function  of  projectile  type  (fan, 
compressor,  and  turbine  blades).  Similarly,  standard  deviation  (SD)  is  also  presented  for  each 
model. 

|E3  Fan  (0.05"  Al)  ■  Fan (0.07"  Al)  DComp.  Blurb. 


FIGURE  30.  PREDICTION  ACCURACY  BY  PREDICTION  MODEL 
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As  seen  in  figure  30  and  below  in  table  2,  all  three  of  the  alternative  equations  provide  improved 
standard  deviation  of  the  mean  prediction  values.  Of  the  three,  the  FAA  energy  equation  (with 
modified  constant)  provided  the  most  accurate  predictions  while  maintaining  the  least  standard 
deviation.  Modification  of  the  constant  was  achieved  by  numerically  curvefitting  the  equation’s 
empirical  constant  to  the  test  data.  Although  the  equation  should  provide  good  results  for  2024 
T3  aluminum,  its  performance  is  unknown  should  predictions  for  other  materials  be  desired. 


TABLE  2.  PREDICTION  ACCURACY  BY  PREDICTION  MODEL 


Fan  Blades 
(0.05"  Al.) 

Fan  Blades 
(0.07"  Al.) 

Compressor 

Blades 

Turbine 

Blades 

All  Test 
Series 

JTCG 

Fragment 

0.89/0.39 

0.80/0.26 

0.90/0.38 

1.02/0.14 

0.88/0.36 

JTCG 

Projectile 

0.89/0.12 

0.81/0.21 

0.94/0.22 

0.96/0.09 

0.87/0.27 

FAA  Energy 
Equation 

1.26/0.28 

1.23/0.22 

1.34/0.29 

1.19/0.22 

1.25/0.26 

FAA  Energy 
Equation  (Mod) 

1.08/0.15 

1.02/0.15 

1.10/0.25 

1.07/0.14 

1.06/0.20 

Of  the  three  alternate  equations,  the  one  which  is  the  most  suitable  for  current  use  is  the 
JTCG/ME  equation  for  projectiles.  Although  the  mean  prediction  accuracy  of  this  equation  is 
similar  to  that  of  the  fragment  equation,  standard  deviation  is  greatly  reduced.  This  is  due  to 
better  accounting  for  petaling  effects  in  the  equations  empirical  constants.  Additionally,  the 
database  of  empirical  constants  for  this  equation  is  more  comprehensive  with  respect  to  common 
aircraft  materials  (including  7075  T6  aluminum  and  2024  T3  aluminum). 

It  should  be  noted  that  the  presented  area  is  still  the  main  dependent  variable  which  accounts  for 
fragment  orientation  within  this  equation.  As  such,  error  induced  by  perimeter  shear  is  still 
present  within  the  projectile  equation. 

6.  SUMMARY. 

The  testing  covered  in  this  report  was  conducted  to  verify  the  effectiveness  of  the  JTCG/ME 
penetration  equations.  In  particular,  the  testing  validated  the  equation’s  ability  in  predicting  the 
residual  velocity  of  representative  jet  engine  blade  fragments. 

Limitations  of  the  equations  were  identified  which  influence  the  accuracy  of  the  predictions.  It 
was  determined  that  prediction  accuracy  was  mainly  dependent  on  the  type  of  perforation  mode 
incurred  during  each  impact. 

Those  penetrations  which  resulted  in  plugging  perforation  of  the  target  were  the  most  accurate,  as 
data  from  plugging  perforations  was  utilized  in  the  development  of  the  penetration  equations. 
However,  prediction  accuracy  was  degraded  when  petaling  was  the  predominant  perforation 
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mode.  In  nearly  all  cases,  petaling  perforation  resulted  in  significant  underprediction  of  the 
fragment’s  residual  velocity. 

The  actual  shear  length  along  the  perimeter  of  the  presented  area  was  suggested  to  be  a  possible 
influential  factor  for  near  edge-on  impacts  of  the  engine  fragments.  Overprediction  of  the 
fragment’s  residual  velocity  may  be  expected  for  these  impacts.  This  is  mainly  due  to  the 
difference  in  fragment  shape  between  the  engine  fragments  and  the  fragments  utilized  in  the 
development  of  the  penetration  equations. 

Other  limitations  included  unknown  degradations  to  the  target  material  which  may  occur  from 
aging.  The  importance  of  presented  area  accuracy  was  stressed  for  impact  conditions  which  are 
close  to  the  ballistic  limit.  Finally,  contribution  of  the  fragment’s  rotational  kinetic  energy  was 
determined  to  be  of  minor  significance. 

7.  CONCLUSIONS  AND  RECOMMENDATIONS. 

From  the  compilation  of  the  results  of  all  tests,  the  JTCG/ME  penetration  equations  performed 
adequately.  However,  the  particular  prediction  limitations  of  the  ballistic  limit  equation  for 
fragments  must  be  noted.  As  discussed,  these  limitations  are  mainly  due  to  petaling  perforation 
and  perimeter  shear  induced  errors. 

Of  the  alternative  ballistic  limit  equations  considered,  the  JTCG/ME  ballistic  limit  equation  for 
projectiles  proved  the  most  suitable  for  current  use.  This  is  due  to  the  equation’s  better 
accounting  for  petaling  induced  errors,  which  effectively  reduced  the  standard  deviation  of  mean 
prediction  accuracy.  Additionally,  the  equation  offers  a  more  comprehensive  empirical  constant 
database  of  common  aircraft  materials. 

The  FAA  Energy  Equation  (with  modified  constant)  proved  the  most  accurate  of  all  the 
alternative  ballistic  limit  equations.  This  was  achieved  through  numerical  curvefitting  of  the 
equation’s  empirical  constant  to  the  test  data  (for  2024  T3  aluminum).  The  equation  also  offers 
the  benefit  of  conservative  overpredictions  of  the  residual  velocity.  This  equation  is 
recommended  for  current  use  with  2024  T3  aluminum  only.  It  is  also  recommended  that  its  use 
be  expanded  as  test  data  (and  thus  empirical  constants)  become  available  for  other  common 
aircraft  materials. 

To  ease  analysis  and  evaluation  of  the  penetration  equations,  all  of  the  tests  were  conducted  with 
normal  impacts  of  the  targets.  It  is  currently  unknown  if  oblique  target  impacts  would  provide 
the  same  results.  It  is  thus  recommended  that  some  future  testing  be  conducted  to  evaluate  these 
equation’s  prediction  performance  in  oblique  impacts. 

A  high  degree  of  unintentional  stringer  impacts  occurred  during  engine  cowling  testing.  This 
suggests  that  aircraft  structural  members  (other  than  skin)  may  also  play  an  integral  part  in  engine 
debris  damage  mitigation.  It  is  recommended  that  these  structures  also  be  included  in  the 
Uncontained  Debris  Model.  It  is  also  recommended  that  some  testing  be  conducted  to  further 
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refine  the  role  these  structures  play  in  damage  mitigation.  The  penetration  equations  may  thus  be 
refined  through  this  testing  to  enable  post-stringer  impact  residual  velocity  predictions. 

A  metallurgical  analysis  should  be  conducted  on  the  cowl  material  to  determine  the  effect 
corrosion  and/or  aging  may  have  had  on  the  cowl  door  material  properties  and  penetration  test 
results.  The  results  of  this  analysis  could  be  reported  in  a  subsequent  report. 
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APPENDIX  A— DERIVATION  OF  THE  JTCG/ME  PENETRATION  EQUATIONS 

The  Ballistic  Limit  Equation 


The  JTCG/ME  Penetration  Equations  include  a  relationship  for  calculating  the  ballistic  limit 
velocity  (V50).  This  relation  may  be  written  for  a  non-oblique  impact  as  follows  (subscript  “n” 
denotes  non-oblique,  or  normal  impact): 


ViOn  —  Q.f 


PMp)  b' 

W  J 


PdApV 
W0  J 


(1) 


Equation  (1)  is  a  non-dimensionalized  form  of  the  following  empirical  relation: 


VjOn  =  10c(tA  J(W)"P  (2) 

Equation  (2)  is  one  of  a  series  of  ballistics  relations  developed  by  Johns  Hopkins  University  for 
the  Ballistics  Research  Laboratory  under  Project  THOR.  It  is  an  empirical  relation  describing  the 
ballistic  limit  velocity  for  warhead  fragments  of  an  unspecified  shape  impacting  a  metallic  target. 


Equation  (2)  includes  three  empirically  derived  constants  a,  p,  and  c  for  specified  target 
materials.  Further  examination  reveals  there  are  only  three  dependent  variables.  These  are  target 
thickness  (t),  weight  of  the  fragment  (W),  and  presented  area  of  the  fragment  to  the  target  plate 
(Ap).  An  interpretation  of  these  variables  may  be  made  such  that  for  a  specified  target  material, 
the  only  factors  affecting  the  ballistic  limit  velocity  are  the  thickness  of  the  plate,  the  mass  of  the 
impactor,  and  the  area  of  the  target  in  which  the  impactor’s  kinetic  energy  is  distributed. 

It  should  be  stressed  that  equation  (2)  requires  the  strict  use  of  specified  units  for  the  variables. 
Thickness  and  presented  area  must  be  in  units  of  inches  and  square  inches,  respectively. 
Fragment  weight  must  be  in  units  of  grains  weight  (7000  grains  =  1  Pound).  The  V50  output  of 
the  equation  is  in  units  of  feet  per  second. 

To  reduce  the  possibility  of  incurring  errors  from  the  input  of  incorrect  units  into  equation  (2), 
equation  (1)  was  developed  in  an  attempt  to  non-dimensionalize  these  inputs.  Equation  (2)  may 
be  rewritten  as: 


(3) 


The  right  two  groupings  of  terms  in  equation  (3)  are  dimensionless.  Two  new  terms  have  been 
introduced  to  achieve  these  dimensionless  groupings,  but  they  have  no  influence  (as  they  cancel) 
on  the  results  of  the  equation.  These  terms  include  fragment  specific  density  (pf)  and  a  specified 
reference  weight  constant  (Wo  =100  grains). 

Comparison  of  equation  (3)  with  equation  (1)  reveals  that  the  empirical  constants  utilized  by  the 
JTCG/ME  equation  are  in  fact 
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br  =  p 
f  =  a-  p 
W, 
P 


cbf  = 


r  \  a 

fj  (Wo)" 


10c 


We  note  that  bf,  and  f  are  indeed  constants,  but  that  Cbf  is  in  fact  a  variable.  As  W0,  a,  p,  and  c 
are  all  constants,  Cbf  is  dependent  only  on  the  fragment’s  specific  density  Pf. 

In  summary,  equation  (2)  may  be  utilized  with  any  type  or  material  of  fragment  for  a  specified 
target  material  (as  defined  by  empirical  constants).  However,  the  use  of  equation  (1)  requires  the 
conversion  of  Cbf  for  fragments  of  different  materials  (or  specific  densities). 

The  Residual  Velocity  Equation 

Derivation  of  the  residual  velocity  equation  is  roughly  outlined  on  pages  83-84  of  the  JTCG/ME 
Penetration  Equations  Handbook  (61  JTCG/ME-77-16).  However,  many  steps  of  this  derivation 
are  omitted  within  the  handbook  due  to  space  constraints.  The  derivation  of  this  equation  is 
presented  here  in  its  entirety  to  illustrate  the  assumptions  on  which  this  equation  is  based.  As 

such,  it  is  hoped  that  users  of  this  equation  will  more  readily  understand  the  limitations  of  the 
equation. 

The  derivation  is  achieved  through  the  application  of  the  conservation  of  energy  and  momentum 
for  ballistic  impact  of  a  blunt  projectile.  The  main  assumptions  utilized  in  this  derivation  are 

Assumption  #1:  Assume  that  a  plugging  type  of  penetration  occurs.  At  penetration,  a  plug  is 
sheared  from  the  target  and  is  propelled  forward  with  the  penetrator  along  the 
penetrators  original  trajectory.  The  velocity  of  the  plug  will  be  at  the  same 
residual  velocity  as  the  penetrator. 

Assumption  #2:  Assume  that  if  penetration  occurs  for  a  given  projectile  at  a  specified  impact 
orientation,  the  energy  dissipated  at  penetration  is  a  constant  regardless  of 
impact  velocity.  A  value  for  this  energy  loss  may  be  determined  as  being  equal 
to  the  kinetic  energy  of  the  fragment  at  its  ballistic  limit  (i.e.,  V50). 

Assumption  #3 :  Assume  that  the  plate  is  much  greater  in  mass  than  the  impacting  projectile. 

Assumption  #4:  It  is  assumed  that  penetration  of  the  plate  occurs  in  two  distinct  phases: 


Phase  1: 


The  penetrator  is  traveling  at  a  specified  impact  velocity.  It  then  impacts  the 
target  plate,  resulting  in  a  shearing  of  the  plate  which  further  results  in  the 
formation  of  a  plug  of  plate  material.  The  bounding  area  of  this  plug  is 
equivalent  to  the  penetrator’s  presented  area  at  the  moment  of  impact.  The 
phase  ends  just  as  this  shearing  is  completed.  At  this  completion,  the 
penetrator  has  been  slowed  to  a  postimpact  velocity  denoted  as  V(bar).  The 
target  plate  has  acquired  some  of  the  penetrator’s  momentum  and  is  now 
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traveling  at  a  postimpact  velocity  of  Vt.  The  plug  has  been  formed  but  has  not 
yet  accelerated. 

Phase  2:  Upon  shearing  of  the  plug  from  the  plate,  the  plug  and  fragment  undergo 

acceleration  and  deceleration  respectively.  The  phase  ends  as  both  the  plug  and 
fragment  achieve  an  equivalent  residual  velocity. 


Figure  A-l.  Derivation  of  the  Residual  Velocity  Equation 


Phase  1. 

During  phase  1,  the  fragment  is  initially  traveling  at  an  impact  velocity  (V).  It  then  impacts  the 
plate  resulting  in  momentum  transfer  to  the  plate. 

By  the  conservation  of  momentum: 

Fragment  Impact  Momentum  =  Fragment  Post-Impact  Momentum  +  Plate  Post-Impact  Momentum 


mV  =  mV  +  mtVt 

or  v,  =— (v-v) 

mtv  y 

At  impact,  some  of  the  fragment’s  kinetic  energy  is  transferred  to  the  plate.  As  such,  the  plate 
acquires  some  kinetic  energy  of  its  own.  Also  some  of  the  fragment’s  kinetic  energy  is  dissipated 
through  deformation  and  shearing  of  the  plate. 

By  the  conservation  of  energy: 

Fragment  Impact  KE  =  Fragment  Post-Impact  KE  +  Plate  Post-Impact  KE  +  Energy  Dissipation 

~  mV2  =  mV2  +  mV2  +  E 
2  2  2  1 
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But  the  plate  is  of  much  greater  mass  than  the  fragment.  Although  the  plate  did  acquire 
momentum  from  the  impact,  its  acquired  kinetic  energy  is  negligible  in  comparison  with  the 
other  terms  in  the  energy  balance.  This  is  mainly  due  to  the  squared  velocity  terms  in  the 
equation.  The  energy  equation  may  be  rewritten  as 

jmV2  =  imV2  +  E0 

Now  consider  the  dissipated  energy  term.  By  assumption  #  2,  for  a  given  fragment  and  impact 
orientation,  the  energy  dissipated  by  shear  and  deformation  is  a  constant  for  all  impacts  above  the 
ballistic  limit.  In  effect,  if  the  target  plate  fails  at  one  velocity,  impact  at  higher  velocities  will 
also  result  in  failure.  However,  the  energy  dissipated  by  these  failures  is  assumed  to  be  constant. 
It  is  further  assumed  that  this  dissipated  energy  is  equivalent  to  the  fragment’s  kinetic  energy  at 
the  ballistic  limit.  As  such,  the  energy  equation  is  further  modified  to  read 

2  mV2  =  ~  mV2  +  ^  mV52 
or 

V=Vv2-V/0 


Phase  2. 

Phase  2  commences  just  after  formation  and  release  of  the  plug  from  the  plate.  The  plug 
undergoes  acceleration  and  deformation  during  this  phase.  Additionally,  the  fragment  undergoes 
deceleration.  The  phase  ends  just  as  the  plug’s  velocity  equals  the  fragments  velocity.  The  plug 
and  fragment  then  continue  onward  at  this  final  residual  velocity. 

By  the  conservation  of  momentum  for  phase  2: 


Fragment  Post-Impact  Momentum  -  Fragment  Residual  Momentum  +  Plate  Residual  Momentum 


mV  =  mV  +  mpVr 

But  from  phase  1,  we  know  v=7v2-y£.  Substituting  this  relation  into  the  conservation  of 
momentum  for  phase  2  yields: 


Vv2-V50 

mn 


But  the  mass  of  the  plug  is  simply  the  volume  of  the  plug  times  its  specific  weight.  Also  the 
volume  of  the  plug  is  equal  to  the  area  of  the  plug  times  its  thickness.  Note  that  the  plug’s  area  is 
also  equal  to  the  presented  area  of  the  fragment  for  normal  impacts.  However,  for  oblique 
impacts,  the  plug’s  area  is  equal  to  the  fragment’s  presented  area  divided  by  the  cosine  of  the 
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obliquity  angle.  Substituting  into  the  above  relation  yields  the  final  form  of  the  residual  velocity 
equation: 


1+ — — 

WcosG 

It  should  be  noted  that  during  the  acceleration  of  the  plug  and  deceleration  of  the  fragment, 
energy  is  dissipated  through  deformation  of  the  plug.  By  the  conservation  of  energy  for  phase  2: 

Fragment  Post-Impact  KE  =  Fragment  Residual  KE  +  Plate  Residual  KE  +  Energy  Dissipation 

imVi=imVr2  +  impVr2  +  Ep 
The  energy  dissipated  through  deformation  of  the  plug  is  thus: 
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Test  Series  2:  Simulated  Blade  Fragment  vs.  0.05"  Aluminum  Target 


APPENDIX  B— COMPILATION  OF  TEST  RESULTS 
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Test  Series  3:  Fan  Blade  Fragment  vs.  0.05"  Aluminum  Target 
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Shot  Number 

Figure  B-5.  Test  Series  3,  Prediction  Accuracy 
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Figure  B-6.  Test  Series  3,  Impact  and  Residual  Velocity  Median  Comparison 
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Test  Series  4:  Fan  Blade  Fragment  vs.  0.071" Aluminum  Target 
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Figure  B-9.  Test  Series  4,  Impact  and  Residual  Velocity  Median  Comparison 
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Test  Series  5  and  7:  Compressor  Blade  Fragment  vs.  0.05 "Aluminum  Target 
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Figure  B-10.  Test  Series  5  and  7,  Impact  and  Residual  Velocity  Comparison 
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Figure  B-ll.  Test  Series  5  and  7,  Prediction  Accuracy 
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Figure  B-12.  Test  Series  5  and  7,  Impact  and  Residual  Velocity  Median  Comparison 
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Test  Series  6:  Turbine  Blade  Fragment  vs.  0.05 " Aluminum  Target 
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Test  Series  8:  Fan  Blade  Fragment  vs.  Aircraft  Engine  Cowling 
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Figure  B-17.  Test  Series  8,  Prediction  Accuracy 
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Figure  B-18.  Test  Series  8,  Impact  and  Residual  Velocity  Median  Comparison 
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APPENDIX  C— TARGET  IMPACT  PHOTOGRAPHS 


Preface 

This  appendix  has  been  included  to  present  the  test  data  in  an  easily  interpreted  form.  The  intent 
is  to  present  data  which  may  be  most  useful  for  follow-on  damage  prediction  efforts,  in  particular 
for  finite  element  analysis. 

The  data  is  tabulated  as  follows: 

Target/Thickness:  The  material  composition  of  the  target  and  its  thickness  in  inches. 

Fragment/Material/Weight:  The  type  of  fragment  utilized  (Fan,  turbine,  compressor,  or 
simulated  fragment)  along  with  its  material  composition  and  weight  in  grams. 

Length/Width/Thickness:  The  length,  width,  and  thickness  of  the  fragment  in  inches.  The 
dimensions  describe  an  equivalent  solid  parallelepiped  (rectangle  of  constant  thickness)  of  the 
same  material  and  weight  as  the  actual  test  fragment.  The  listed  dimensions  for  thickness  are 
average,  as  most  of  the  fragments  utilized  were  in  the  shape  of  airfoils  with  varying  thickness 
along  the  length  of  the  fragment. 

Pitch/Yaw/Roll:  The  pre-impact  orientation  of  the  fragment  in  terms  of  sequential  axis 

rotations  of  pitch,  yaw,  and  roll  (degrees).  For  amplified  discussion,  please  refer  to  the 
orientation  discussion  of  the  main  report. 

Presented  Area:  The  pre-impact  presented  area  of  the  fragment  as  projected  along  the 

fragment’s  trajectory  onto  the  target. 

Velocity  Comparison  (Impact/Residual/Predicted/Accuracy):  The  pre-impact  and  post¬ 
penetration  residual  velocity  (in  feet/second)  of  the  fragment  as  measured  by  high-speed  film 
footage.  For  comparison  purposes,  the  residual  velocity  as  predicted  by  use  of  the  JTCG/ME 
penetration  equations  (for  fragments)  has  also  been  included,  along  with  this  value’s  prediction 
accuracy  (in  comparison  with  the  actual  residual  velocity). 

* 

Additional  data  is  presented  via  included  figures  for  each  shot.  These  include: 

Pre-Impact  Orientation:  This  figure  is  a  computer  rendered  illustration  of  the  fragment  and  its 
orientation  at  the  moment  just  prior  to  impact  with  the  target.  The  fragment  is  centered  within  a 
three-dimensional  axis  describing  the  pre-rotation  reference  pitch  (to  the  left),  yaw  (up),  and  roll 
(into  the  target)  axes. 

Three  shadow  projecting  light  sources  have  been  included  in  the  illustration.  These  include  a 
light  source  aimed  along  the  trajectory  of  the  fragment  and  projecting  a  shadow  of  the  fragment’s 
presented  area  onto  the  target.  Two  additional  light  sources  are  included  which  project 


C-l 


horizontal  and  vertical  plane  shadows  of  the  fragment,  simulating  pre-impact  views  of  the 
fragment  as  seen  on  the  high-speed  camera  footage. 

A  1/2"  square  checkerboard  pattern  has  been  utilized  on  the  target,  floor,  and  wall  planes.  This 
pattern  has  been  included  as  a  visual  dimensional  reference.  In  addition,  the  target  perimeter  has 
been  outlined  with  a  darkened  checkerboard  pattern  simulating  the  7.5"  x  9.5"  photographic 
gridboard  described  below  and  included  within  the  actual  damage  photographs  of  each  target. 

Target  Damage  (Frontal  and  Rear  Views):  These  figures  are  actual  photographs  of  the  post¬ 
impact  damage  incurred  by  each  target.  The  damage  has  been  correctly  oriented  and  centered 
within  a  7.5"  x  9.5"  photographic  gridboard  of  1/2"  checkered  squares.  Note  that  each 
photograph  has  been  cropped  along  the  edge  of  the  gridboard.  This  was  done  to  allow 
enlargement  of  each  photo  to  emphasize  the  characteristics  of  the  impact  damage.  As  such  the 
entire  area  of  the  target  (2'  x  2')  is  not  visible  within  each  picture.  However,  in  nearly  all  test 
shots,  the  damage  noted  outside  of  the  gridboard  was  negligible. 
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Figure  C-110:  Target  Damage  Figure  C-lll:  Target  Damage 

Frontal  View  Rear  View 
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Background 

As  part  of  the  ongoing  Uncontained  Engine  Debris  Mitigation  Program,  NAWCWPNS 
has  developed  a  model  to  calculate  uncontained  engine  debris  exit  velocities  and  residual 
energy  levels.  The  Uncontained  Debris  Model  (UDM)  is  based  on  the  Joint  Technical 
Coordinating  Group  for  Munitions  Effectiveness  (JTCG/ME)  penetration  equations 
which  have  been  derived  from  50  caliber  projectile  testing.  Engine  debris  characteristics 
are  significantly  different  from  the  50  cal.  projectile  in  that  the  velocity  range  is  lower, 
debris  fragment  size  is  larger,  and  the  mass  is  heavier. 

The  UDM  focuses  on  blade  fragment  debris  as  this  type  of  debris  is  most  likely  to  result 
from  an  uncontained  engine  event.  According  to  extensive  investigation,  this  is  the  most 
likely  particle  which  would  result  in  damage  to  aircraft.  The  debris  model  identifies  the 
debris  types  by  component  (fan  blade,  turbine  blade, ...)  and  defines  the  most  likely  size, 
mass  and  velocity  for  each  debris  fragment. 

Debris  fragment  penetration  orientation  is  defined  for  each  impact.  The  fragment’s 
impact  orientation  angle  is  defined  as  the  angle  between  the  fragment’s  surface  area  and 
the  perpendicular  of  the  fragment’s  direction  of  travel.  Obliquity  angle  is  defined  as  the 
angle  between  the  fragment’s  direction  of  travel  and  the  perpendicular  of  the  impact 
target  surface  area.  An  illustration  defining  impact  orientation  and  obliquity  is  given  in 
figure  1 .  For  further  clarification,  examples  of  various  orientation  and  obliquity  angle 
impacts  are  shown  in  figure  2. 


Figure  1:  Impact  Orientation  and  Obliquity  Angles 
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Impact  Surface 


\  \  \ 

30  degree  Obliquity  Angle  30  degree  Obliquity  Angle  30  degree  Obliquity  Angle 
90  degree  Impact  Orientation  30  degree  Impact  Orientation  0  degree  Impact  Orientation 


Figure  2:  Examples  of  Obliquity  and  Impact  Orientation 

For  initial  impacts  it  is  assumed  (based  on  investigative  data)  that  as  the  debris  is 
penetrating  case  structure  or  inner  cowl  barrel  walls,  the  impact  angle  is  large  (a  90 
degree  angle  defines  an  end  on  penetration).  Penetrations  through  additional  structures 
are  assumed  to  take  place  at  a  45  degree  impact  orientation  (this  replicates  the  average 
impact  orientation  of  a  tumbling  fragment). 


The  JTCG  penetration  equations  utilized  by  the  model  are  listed  as  follows: 


1. 


1.0  +  - 


PApt 


W  cos# 


V50  "  Cbf 


w 


w 


Variables:  V 

V50 

time. 

V, 

Pf 

P 

t 

4 

e 

w 

W0 


-  debris  initial  velocity 

=  debris  ballistic  limit  velocity.  Penetration  occurs  50%  of  the 

=  debris  residual  velocity 
=  debris  specific  weight  (debris  weight/volume) 

=  plate  specific  weight  (plate  weight/volume) 

=  plate  thickness 

=  debris  presented  area  along  direction  of  travel 

=  debris  obliquity  angle  at  impact 
=  debris  weight 
=  100  grains 


Empirical  Constants:  Cbf,bf,h,f 
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Of  particular  concern  (and  the  reason  for  this  test  series)  are  the  empirical  constants 
utilized  within  equation  2.  Published  values  for  these  constants  may  not  be  appropriate 
for  use  in  the  model.  The  projectiles  which  were  utilized  in  the  development  of  these 
constants  may  have  characteristics  which  are  dissimilar  from  engine  blade  fragments.  As 
such,  prediction  inaccuracies  may  result  from  their  use.  This  test  series  will  verify  the 
accuracy  and  appropriateness  of  these  constants,  and  if  necessary  determine  new  values 
for  these  constants. 


Objective 

The  test  objective  is  to  collect  fragment  penetration  data  to  validate  and/or  refine  the 
penetration  equations  used  to  develop  the  Uncontained  Debris  Model.  Testing  will  be 
conducted  using  a  fragment  simulator  and  engine  debris  fragments  against  various 
thickness’  of  aluminum  plates  and  possibly  aircraft  skin  structures. 

Approach 

Test  Overview 

Fragment  penetration  data  will  be  collected  utilizing  the  setup  displayed  in  figure  3.  An 
airgun  will  be  utilized  to  propel  a  sabot  and  fragment  to  the  desired  velocity.  A  sabot 
stopper  will  be  affixed  to  the  end  of  the  gun  barrel  to  strip  the  sabots  away  from  the 
accelerated  fragment.  Blast  deflectors  will  be  placed  in  front  of  the  gun  muzzle,  and 
beside  the  high  speed  cameras  to  reduce  the  potential  of  blast  damage  to  the  target  and 
cameras. 

Impact  projectile  velocity  will  be  computed  based  on  input  from  a  breakwire  affixed  to 
the  gun  muzzle,  and  breakpaper  affixed  to  the  target  impact  point.  Residual  projectile 
velocity  will  be  computed  based  on  input  from  two  breakpaper  stands  located  at  the  rear 
of  the  target.  Gridboards  will  be  placed  near  the  fragment-target  impact  point  to  allow 
additional  velocity  calculations  to  be  made  from  high  speed  camera  footage.  A  mirror 
will  be  placed  above  the  target  to  assist  determination  of  fragment  orientation  at  impact. 
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Top  View  High  Speed 

Video 

Entrance/Exit 


Figure  3:  Test  Setup 
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The  airgun  will  propel  the  fragments  utilizing  a  sabot  as  displayed  in  figure  4.  The 
desired  fragment  orientation  angle  will  be  controlled  through  a  combination  of  fragment 
tilt  angles  within  the  sabot  (see  figure  5),  and  the  overall  distance  from  the  gun  muzzle  to 
the  target. 


Insert 


Side  Front 

Figure  4:  Saboted  Fragment  (90  deg.  impact  orientation) 


Insert 

Side  Front 

Figure  5:  Sabot  Impact  Orientation  Angle  Control,  (60  degrees  shown) 
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The  fragment  will  be  released  from  the  sabot  through  the  use  of  a  sabot  catcher,  as 
displayed  in  figure  6.  The  tapered  reduction  in  internal  diameter  of  the  catcher  enables 
minimized  impulsive  forces  to  be  seen  by  the  sabot.  The  intent  is  to  reduce  the 
possibility  of  breakup,  and  travel  of  sabot  fragments  to  the  target.  However,  the  sabots 
are  not  reusable  due  to  severe  scoring  damage  from  the  catcher. 


Blast 

Exhauster 


Figure  6:  Sabot  Catcher 


The  sabot  catcher  will  be  affixed  a  distance  from  the  air  gun  muzzle  by  means  of  a  blast 
exhauster.  The  lodging  of  the  sabot  within  the  catcher  will  reduce  blast  exhaust  through 
the  catcher  exit  (and  thus  blast  seen  by  the  target).  The  resulting  pressure  build  up  behind 
the  sabot  will  be  released  through  the  exhauster.  The  exhauster  will  allow  much  of  the 
blast  to  dissipate  perpendicular  to  the  muzzle,  thus  minimizing  blast  at  the  target. 


D-8 


6 


Test  Sequence 

The  test  will  be  composed  of  eight  separate  test  series  as  outlined  in  Table  1 .  A  more 
detailed  test  matrix  is  outlined  in  the  Test  Procedure  section  of  this  test  plan  as  Table  4. 


Table  1:  Test  Sequence 


Test 

Series 

Fragment 

(Impact  Orientation) 

Approximate 
Fragment  Weight 

Target 

Number  of 
Tests 

1 

Sabot  Checkout 

0.25  lb. 

0.05”  A1 

3 

2 

Simulated  Frag.  (90) 

0.055  lb. 

0.05”  A1 

8 

3 

Fan  Blade  (90) 

0.4  lb. 

0.05”  A1 

8 

4 

Fan  Blade  (90) 

0.4  lb. 

0.07”  A1 

8 

5 

Compressor  Blade  (45) 

0.251b. 

0.05”  A1 

8 

6 

Turbine  Blade  (60) 

0.25  lb. 

0.05”  A1 

8 

7 

Turbine  Blade  (30) 

0.25  lb. 

o 

o 

Uj 

> 

8 

8 

Fan  Blade  (90) 

0.4  lb. 

Engine  Cowl 

6 

Test  series  #1  will  verify  the  effectiveness  of  the  sabots  in  achieving  both  end  on  and 
tumbled  impacts  at  the  target.  Proper  operation  of  all  test  systems  will  also  be  verified 
during  these  tests.  Additionally,  the  accuracy  of  fragment  velocity  calculations  will  be 
compared  based  upon  break  wire-break  paper  computation  and  high  speed  camera  visual 
calculations. 

Test  Series  #2  will  fire  simulated  blade  fragments  similar  to  those  utilized  during  the  SRI 
International  Advanced  Armor  Technology  impact  experiments.  These  experiments  were 
conducted  in  the  summer  of  1997  at  the  SRI  International  gas  gun  facility.  The  SRI 
experiments  evaluated  the  impact  resistance  of  current  fuselage  wall  materials  (including 
2024-T3  Aluminum),  along  with  several  advanced  candidate  wall  materials.  Test  Series 
#2  will  provide  baseline  impact  resistance  data  for  2024-T3  Aluminum.  This  baseline 
data  (along  with  the  Uncontained  Debris  Model  predictions)  will  be  compared  with  the 
accumulated  data  from  the  SRI  experiments. 

Test  series  #  3  through  7  will  examine  actual  engine  blade  fragment  impacts  at  various 
impact  orientation  and  obliquity  angles.  Characteristics  of  these  fragments  (dimensions, 
material,  weight,  and  impact  velocity/orientation)  will  be  those  which  are  most  likely  to 
be  experienced  during  an  engine  burst. 

Test  series  #  8  is  currently  planned  as  the  final  test  series,  and  will  examine  the  reaction 
of  engine  blade  impacts  against  an  actual  engine  cowling.  Comparison  will  be  made  with 
the  data  accumulated  from  the  previously  conducted  test  series.  The  intent  is  to 
determine  consistency  of  the  data  between  the  use  of  new  sheet  metal  targets,  and  actual 
engine  cowling  material.  Depending  on  cost,  range  and  material  availability,  an 
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additional  test  series  may  be  included  which  will  examine  fragment  impacts  against 
Kevlar  target  plates.  Details  on  the  Kevlar  plate  test  series  will  be  supplied  by  SRI 
International. 

Fragment  impact  velocities  utilized  in  each  test  series  will  be  based  on  the  predicted 
fragment  V50  values  for  the  test  conditions.  Due  to  the  uncertainty  of  the  actual  V50  value 
for  each  test  condition,  impact  velocities  near  the  predicted  V50  value  may  result  in  non¬ 
penetration  of  the  target  (and  thus  loss  of  desired  residual  velocity  data).  As  such,  the 
initial  test  velocity  utilized  in  each  series  will  be  the  fragment’s  predicted  V50  velocity 
plus  500  ft/sec.  Follow  on  test  impact  velocities  will  incrementally  decrease  from  this 
initial  value  by  V50  plus  300, 100,  75,  50,  25,  and  15  ft/sec.  The  final  test  series  impact 
velocity  will  be  the  predicted  value  of  V50. 

It  should  be  noted  that  there  are  not  enough  shots  available  in  each  test  series  to 
determine  a  precise  value  of  V50  for  each  test  condition.  Should  non-penetration  of  the 
target  occur  prior  to  reaching  the  predicted  value  of  V50,  that  impact  velocity  which 
resulted  in  non-penetration  will  be  assumed  as  an  approximate  value  for  V50.  Those 
remaining  shots  in  the  test  series  will  be  incremented  from  this  approximate  value. 

Test  series  #  2,  3,  and  4  will  examine  a  90  0  (edge  on)  fragment  impact  orientation.  The 
predicted  V50  value  for  this  type  of  condition  can  be  quite  low  (<  200  ft/sec)  due  to  the 
small  fragment  presented  area  (assuming  each  fragment  is  roughly  a  flat  plate).  Very  low 
impact  velocities  (<  200  ft/sec)  may  prove  difficult  to  accurately  regulate  with  the  airgun. 
These  velocities  will  be  attempted,  but  the  actual  impact  velocity  may  vary  from  those 
utilized  in  the  pre-test  predictions. 

Many  of  the  tests  will  utilize  fragments  cut  from  actual  engine  blades.  As  shown  in 
figure  7,  these  fragments  may  have  a  slight  aerodynamic  twist  in  their  dimensions.  As 
compared  to  a  flat  plate,  a  twist  will  result  in  an  increased  actual  edge-on  presented  area, 
and  thus  increased  V50  values  for  edge-on  impacts.  The  fragment’s  actual  edge-on 
presented  area  will  be  determined  after  fragment  fabrication  (and  prior  to  testing)  to  allow 
refinement  of  the  pre-test  predictions. 

m 


Edge-On  View  Edge-On  View 

Flat  Plate  With  Twist 

Figure  7:  Example  of  Fragment  Edge-On  Views 
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For  impact  velocities  which  result  in  penetration  of  the  target,  the  fragment’s  residual 
velocity  will  be  measured.  Data  will  be  collected  by  computing  the  fragment’s  impact 
velocity  via  a  breakwire  in  the  gun  barrel  and  break  paper  on  the  target.  High  speed  film 
will  be  used  to  verify  the  impact  velocity  and  calculate  the  residual  velocity.  High  speed 
film  will  also  be  used  to  determine  the  fragment’s  impact  orientation.  A  sample  chart 
illustrating  the  type  of  data  to  be  expected  from  this  test  series  is  illustrated  in  figure  8. 

Sample  Data  Chart 

Residual  Velocity  (ft/sec) 

800 

700 
600 
500 
400 
300 
200 
100 
0 

0  200  400  600  800 

Impact  Velocity  (ft/sec) 

Figure  8:  Sample  Velocity  Data  Chart 
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Pretest  Predictions 

The  JTCG/ME  penetration  equations  prediction  model  was  utilized  to  calculate  predicted 
residual  velocities  for  each  test  series.  V50  values  were  predicted  based  on  each  test  series 
parameters.  Predicted  residual  velocities  were  then  tabulated  based  on  impact  velocity 
increments  of  V50  +  15, 25,  50,  75,  100,  300,  and  500  fit/sec.  A  prediction  plot  of  each 
test  series  residual  velocities  is  illustrated  in  figure  9. 


Test  Series  6 

Test  Series  2  Turbine  Blade  (60) 

Sim.  Fragment  (90)  m=0.25  lbs 


Figure  9:  Pre-Test  Predictions,  Initial  Impact  vs.  Residual  Velocity 

The  fragment’s  mass  was  then  utilized  in  conjunction  with  impact  and  residual  velocities 
to  calculate  impact  Kinetic  Energy  versus  Residual  Kinetic  Energy.  Figure  10  displays 
this  comparison. 


Figure  10:  Pre-Test  Prediction,  Initial  vs.  Residual  Kinetic  Energy 
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An  alternative  comparison  was  made  of  the  ratio  of  the  fragment’s  predicted  residual 
kinetic  energy  to  impact  kinetic  energy.  This  comparison  was  plotted  against  impact 
velocity,  and  is  displayed  in  figure  1 1 . 


Figure  11:  Pre-Test  Prediction,  Impact  Velocity  vs.  Kinetic  Energy  Ratio 
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Pre-Test  Preparation  and  Test  Set-Up 

Fabrication  Requirements 

Target  Plates  and  Support  Stand 

Approximately  forty  three  0.05”,  and  eight  0.07”  thick  aluminum  (2024-T3)  target  plates 
will  be  required  for  the  test  series.  The  plates  should  be  cut  as  2  foot  square  panels,  and 
should  include  drilled  holes  for  connection  to  the  test  stand  (see  figure  12). 


Figure  12:  2x2  ft.  Target  Plate 

For  test  series  #  8,  the  required  eight  target  plates  will  be  constructed  of  actual  engine 
cowling  material  from  WSL  aircraft.  An  additional  test  series  (#  9)  may  be  added  which 
would  require  the  preparation  of  Kevlar  targets  utilizing  Kevlar  material  provided  by  the 
WSL.  Test  preparation  details  of  the  Kevlar  targets  will  be  supplied  by  SRI. 

i 

A  target  support  test  stand  will  be  constructed  of  angle  iron  as  displayed  in  figure  13.  The 
stand  will  include  drilled  holes  for  bolting  the  target  plate  to  the  stand.  Four  steel  support 
strips  with  matching  drilled  holes  will  also  be  fabricated  to  provide  rigid  attachment  of 
the  plates  to  the  test  stand. 


Welded 
Angle  Iron 
Support  Frame 


Figure  13:  Plate  Support  Stand 
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Residual  Velocity  Break  Paper  Stands 


Two  residual  velocity  break  paper  support  stands  will  be  required  for  the  test  series. 
Fabrication  will  be  similar  to  that  displayed  in  figure  14. 


2'  X  2' 
Plywood 
Frame 


2"  X  2"  Pine 
Stand 


Figure  14:  Residual  Velocity  Break  Paper  Stands 


Projectile  Fragments 

Test  series  #2  will  require  the  fabrication  of  8  simulated  fan  blade  projectiles.  They  will 
be  fabricated  from  titanium  plate.  Fragment  dimensions  will  be  1  3/8”  x  1”  x  1/4”.  Each 
fragment  will  weigh  approximately  0.055  lbs. 

Test  series  #  3, 4,  and  8  will  utilize  actual  engine  fan  blade  fragments.  They  will  be  cut  to 
dimensions  of  5  1/2”  length  and  3”  width.  Each  fragment  will  weigh  approximately  0.4 
lbs.  A  total  of  22  fan  blade  fragments  will  be  required. 


Test  series  #  5  will  require  the  cutting  and  preparation  of  eight  compressor  blade 
fragments.  Desired  length  and  width  of  the  fragments  are  4”  x  2”.  Each  fragment  will 
weigh  approximately  0.25  lbs. 


Test  series  6  &  7  will  require  preparation  of  16  turbine  blade  fragments,  cut  to 
dimensions  of  4”  x  2”.  Each  fragment  will  weigh  approximately  0.25  lbs. 

A  summary  of  the  required  fragment  projectiles  is  listed  in  Table  2.  Figure  15  illustrates 
how  the  fragments  will  be  cut  from  available  engine  fan,  compressor,  and  turbine  blades. 
Drawings  of  the  represented  fragments  are  displayed  in  Figures  16, 17,  and  18. 
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Table  2:  Projectile  Requirements 


Material 

Dimensions 

L  x  W  x  thickness 

Weight 

Titanium 

1.4”  x  l”x  0.25” 

0.055  lb. 

Fan  Blades 

3”  x  5.5”  x  blade  thickness 

0.4  lb. 

Compressor  Blades 

4”  x  2”  x  blade  thickness 

0.25  lb. 

Turbine  Blades 

4”  x  2”  x  blade  thickness 

0.25  lb. 

#  Fragments 
Required 

8 

22 

8 

16 


I  3-  | 

Fan  Blade 


8"  (Approx.)  — ►If 


K  4  i 

Turbine/ Compressor  Blade 


Figure  15:  Fragment  Preparation 
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1/4 " 

|  |  1  3/8 "  — *i 


FRONT  SIDE 


Titanium,  0.055  lbs 

Figure  16:  Test  Series  2,  Simulated  Fragment 


0.15  " 


FRONT  SIDE 


Titanium,  0.25  lbs 

Figure  17:  Test  Series  5,  6,  &  7,  Compressor  /  Turbine  Blade  Fragment 
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5  1/2 


Figure  18:  Test  Series  3, 4,  &  8,  Fan  Blade  Fragment 
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Sabots  and  Styrofoam  Supports 

Each  of  the  test  series  will  require  sabots  to  propel  the  fragment  projectiles.  Their 
diameter  and  length  should  be  sufficient  to  provide  a  snug  fit  within  the  airgun  barrel, 
while  providing  support  to  the  projectile  fragment.  The  sabots  will  be  fabricated  from 
plastic  tubing  endcaps.  A  total  of  seventy  sabot’s  will  be  required  for  this  testing. 


Styrofoam  insert  supports  will  be  fabricated  and  installed  within  each  sabot  as  previously 
shown  in  figures  4  and  5.  Requirements  for  these  supports  are  listed  in  Table  3. 

Table  3:  Sabot  Styrofoam  Support  Requirements 


Test  Series  Fragment  Width  x  Length 

1  2”  x  4” 

2  1  “  x  1  3/8  “ 

3  3  “  x  5  1/2  “ 

4  3  “x  5  1/2“ 

5  2  ”  x  4  ” 

6  2  ”  x  4  ” 

7  2  ”  x  4  ” 

8  2  ”  x  4  ” 

Sabot  Catcher  and  Blast  Exhauster 


Support  Angle  #  Supports  Required 


45°  3 

90°  (end  on)  8 

90 0  (end  on)  8 

90°  (end  on)  8 

45°  8 

60°  8 

30°  8 

90  0  (end  on)  8 


A  sabot  catcher  /  blast  exhauster  (see  figure  6)  will  be  fabricated  and  installed  to  the 
airgun.  The  catcher  should  enable  easy  removal  of  the  sabot  for  reuse.  Design  and 
materials  utilized  in  this  fabrication  will  be  at  the  discretion  of  the  WSL  Engineer. 


Ordnance  Requirements 

Due  to  the  relatively  large  projectile  size,  testing  will  require  utilization  of  the  Missile 
Impact  Kinetic  Energy  Simulator  (MIKES)  air  gun.  A  high  pressure  nitrogen  gas  source 
will  also  be  required  capable  of  propelling  one  pound  blade  fragments  through  the 
MIKES  at  velocities  up  to  1000  fps. 

Placement  of  the  gun  in  front  of  the  target  will  be  such  that  range  from  the  gun  muzzle  to 
the  target  should  be  within  10  feet.  The  gun’s  shotline  will  be  perpendicular  to  the  target 
plate’s  face  to  achieve  an  obliquity  angle  of  zero  degrees. 
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Target  Setup 


Details  of  the  target  Setup  are  displayed  in  figure  19.  Striped  gridboards  will  be  placed 
in  the  background,  and  below  the  target.  This  will  allow  calculation  of  the  fragment’s  pre 
and  post  impact  velocities  through  examination  of  high  speed  camera  footage.  As 
additional  sources  of  velocity  data,  a  break  wire  will  be  installed  at  the  gun  muzzle,  along 
with  break  papers  affixed  at  the  desired  impact  point  and  at  two  stations  rear  of  the  target 
plate.  A  mirror  will  be  placed  at  a  45  degree  angle  above  the  target  to  enable 
determination  of  the  fragment’s  orientation  at  impact. 


Striped 

Gridboards 


Mirror 


Side  View 


Target 

Plate 


Impact 
Break 
Paper  #1 

Support 
Stand 


Residual  Velocity 
Breakpapers  #2 
&  #3  with 
Support  Stands 


Front  View  Top  View  (Mirror  Removed) 


Figure  19:  Target  Setup 
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Data  Requirements 

Primary  data  required  for  these  tests  are  the  fragment’s  target  impact  and  residual 
velocities,  along  with  fragment  orientation  at  impact.  The  velocity  data  will  be  computed 
through  the  use  of  breakwires  and  breakpaper.  As  secondary  backup,  these  velocities  will 
also  be  calculated  based  on  fragment  position/time  data  from  the  high  speed  camera 
footage.  The  impact  orientation  of  the  fragment  will  be  determined  by  high  speed 
photographic  record  of  the  impact  in  two  planes  (mirrored  top  view  and  direct  side  view). 

As  a  cost  saving  effort,  initial  tests  will  utilize  one  high  speed  film  and  one  high  speed 
video  camera.  The  desired  field  of  view  for  both  cameras  (film  and  video)  is  shown  in 
figure  20  (a).  If  the  photographic  resolution  proves  unsatisfactory,  a  second  high-speed 
film  camera  will  be  utilized  with  field  of  views  as  shown  in  figure  20  (b).  Desired  high¬ 
speed  camera  frame  rates  will  be  2,000  frames/sec  for  all  tests.  The  intent  is  to  minimize 
any  confusion  or  error  which  may  be  incurred  by  changing  frame  rates  between  the 
seventy-five  planned  shots.  This  setting  should  also  provide  the  desired  minimum  of  4 
frames  of  fragment  travel  at  the  highest  fragment  velocities. 


(a)  Single  Camera  View 


Camera  #  1 
Field  of  View 
3'  x  3‘ 


(b)  Dual  Camera  Views 


Figure  20:  Desired  High-Speed  Camera  Views 
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Instrumentation 

Instrumentation  requirements  are  described  as  follows: 

•  Breakwire  at  gun  muzzle 

•  Breakpaper  at  target  plate  surface 

•  Breakpaper  at  residual  velocity  support  stand  #1 

•  Breakpaper  at  residual  velocity  support  stand  #2 

•  High  speed  camera  #1 

•  High  speed  camera  #2  (if  required) 

•  High  speed  video 

•  Range  safety  video 

•  Air  gun  controlling  and  monitor 

Test  Constraints 

Environmental 

All  testing  will  be  accomplished  when  adequate  light  is  present  for  the  high  speed 
cameras.  The  testing  will  be  stopped  if: 

•  There  is  wind  over  25  knots. 

•  There  is  rain  or  snow  in  the  test  area. 


D-22 


20 


Safety  and  Security 

Safety  equipment  must  be  operational  during  the  test:  This  includes  data  recorders  and 
range  safety  video  cameras.  Personnel  should  wear  coveralls  on  the  test  pad. 

Tests  will  be  conducted  following  standard  operating  procedures  prepared  by  the  WSL 
and  approved  by  the  NAWC  Safety  Office.  At  no  time  during  the  tests  will  project 
requirements  overrule  the  safety  of  personnel  or  present  hazards  to  the  facility.  Any 
situation  which  the  firing  officer  or  project  engineer  feels  threatens  the  safety  of 
personnel,  facility,  or  test  items  will  be  cause  for  a  delay  in  testing.  The  Weapons 
Survivability  Laboratory  (WSL),  Naval  Air  Warfare  Center  (NAWCWPNS),  will 
consolidate  the  technical  requirements  of  the  test  plan  into  a  NAWCWPNS-approved 
operating  procedure.  This  procedure  will  be  used  for  all  test  runs  in  the  test  program. 

Logistical 

Failure  of  critical  test  or  safety  equipment  will  stop  the  test  until  repair  is  done. 

Test  Procedure 

The  detailed  test  matrix  is  listed  in  Table  4.  The  test  logic  sequence  is  illustrated  in 
figure  21 .  The  following  is  a  sequential  outline  of  the  test  procedures: 

1 .  Install/repair  the  target  plate . 

2.  Install  new  break  papers  and  break  wires. 

3.  Realign/sight  air  gun  muzzle  to  hit  target  while  achieving  desired  impact  orientation. 

4.  Sight  high  speed  cameras  and  video.  Verify  sufficient  lighting  for  coverage. 

5.  Load  air  gun  barrel  with  sabot  and  projectile. 

6.  Prepare  air  gun  for  firing. 

7.  Commence  data  recording 

8.  Fire  the  projectile. 

9.  Stop  data  recording 
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Table  4:  Test  Series  Matrix 
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Schedule  and  Documentation 

The  test  is  scheduled  for  January  1998  at  the  Weapons  Survivability  Laboratory.  A  test 
report  will  be  prepared  which  will  include  the  test  objectives,  approach,  test  set-up,  and 
results.  The  report  will  be  submitted  for  review  within  two  months  after  test  completion. 


* 
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